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\0 ABSTRACT 

L^' We report results of the first global three-dimensional (3D) magnetohydrodynamic 

Ph (MHD) simulations of the waves excited in an accretion disc by a rotating star with a 

^ dipole magnetic field misaligned from the star's rotation axis (which is aligned with 

^ the disc axis). The main results are the following: (1) If the magnetosphere of the 

star corotates approximately with the inner disc, then we observe a strong one-armed 

J" bending wave (a warp). This warp corotates with the star and has a maximum am- 



plitude between corotation radius and the radius of the vertical resonance. The disc's 
center of mass can deviate from the equatorial plane up to the distance of ^ O.lr. 
^ However, the effective height of the warp can be larger, ^ 0.3r due to the fi- 

— ' nite thickness of the disc. Stars with a range of misalignment angles excite warps. 

^ However, the amplitude of the warps is larger for misalignment angles between 15 

and 60 degrees. The location and amplitude of the warp does not depend on viscos- 
ity, at least for relatively small values of the standard alpha-parameter, up to 0.08. 
\^ (2) If the magnetosphere rotates slower, than the inner disc, then a bending wave is 

^ excited at the disc-magnetosphere boundary, but does not form a large-scale warp. 

^ Instead, persistent, high-frequency oscillations become strong at the inner region of 

the disc. These are (a) trapped density waves which form inside the radius where the 
disc angular velocity has a maximum, and (b) inner bending waves which appear in 
^ the case of accretion through magnetic Raleigh-Taylor instability. These two types of 

waves are connected with the inner disc and their frequencies will vary with accretion 
rate. Bending oscillations at lower frequencies are also excited including global oscil- 
lations of the disc. In cases where the simulation region is small, slowly-precessing 
warp forms with the maximum amplitude at the vertical resonance. The present simu- 
lations are applicable to young stars, cataclysmic variables, and accreting millisecond 
d pulsars. A large-amplitude warp of an accretion disc can periodically obscure the light 

from the star. Different types of waves can be responsible for both the high and low- 
frequency quasi-periodic oscillations (QPOs) observed in different types of stars. Inner 
disc waves can also leave an imprint on frequencies observed in moving hot spots on 
the surface of the star. 

Key words: accretion, dipole — plasmas — magnetic fields — stars. 
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1 INTRODUCTION 

Different types of accreting stars have dynamically important 



magnetic fields, such as young T Tauri stars (e.g., Bouvier et 
aL][2007ai), accreting millisecond pulsars (e.g., van der Klis 



2000 



2004). 



2006| ), and accreting white dwarfs (e.g., ,Warner|1995 
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If a rotating star has a tilted dipole (or other non- 
axisymmetric) magnetic field then it applies a periodic force 
on the inner part of the accretion disc and different types of 
waves are excited and propagate to larger distances (e.g., |Lai| 
|1999| ) . The linear theory of waves in discs has been develop- 
ing extensively over many years (see, e.g., book by |Kato et al.| 
[r998 and references therein). If an external force is applied 
to the disc from the tilted dipole of a rotating star (or from a 
secondary star or planet), then this force can generate density 
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Figure 1. An example of the magnetospheric accretion and a warp 
driven by the rotating and tilted at ^ = 30° dipole magnetosphere 
with p. = 1.5 (model FW/il.5). The background shows one density 
level ip = 0.07). The lines are sample magnetic field lines of the 
closed magnetosphere. The vectors p and show directions of the 
magnetic and rotational moments of the star. 



waves in the disc (e.g., Goldreich||1978|) and vertical bend- 
ing waves (e.g., |Lubow|1981[|Lai|1999| [Ogilvie 2002). Theo- 
retical investigations have led to an understanding that both 
types of waves may be strongly enhanced or damped at radii 
corresponding to resonances (e.g., Kato et ar]|1998| ). There 
are two types of resonances: horizontal resonance (Lindblad 
resonance) and vertical resonance. Another important reso- 
nance is associated with corotation radius rev where the an- 
gular velocity of the disc matches the angular velocity of the 
external force. 

lAlyl (I1980D and |Lipunov| ( [1980| ) investigated the tilting 
of the inner disc due to the magnetic force applied by the ro- 
tating dipole magnetosphere of the star where both the mag- 
netic and rotational axes of the star are misaligned relative to 
the rotational axis of the outer disc. They concluded that the 
magnetic force acting on the inner disc causes it to strongly 
depart from the original plane of the disc. |Lai| ( [1999| ) inves- 
tigated this case and concluded that the magnetic force also 
drives a warp at the inner disc which slowly precesses around 
the rotational axis of the disc. 

[Terquem & Papaloizou] ( [2000l ) investigated the formation 
of a disc warp in the case where the rotational axis of the 
star is aligned with that of the disc, and the inner disc coro- 
tates with the star, taking into account an a— type viscosity 
in the disc ( Shak ura & Sunyaev|l973| . They concluded that 
for a relatively small viscosity, a — 0.01 — 0.08, a warp is 
excited near the magnetosphere and corotates with the star. 
However, it is suppressed by the viscosity at larger distances 
from the star. They also concluded that the height of the warp 
can reach up to ^ 10% of the radial distance. One of the 
main motivations for studying inner disc warps comes from 
observations of young T Tauri stars (CTTS) in which the light 
curves show regular dips that can be interpreted as obscura- 
tion events by a warped disc (e.g., Bouvier et al.|1999[|200"3l 
|2007b| also [Carpenter et al.|200T] ). 

Another important impetus to study waves in the disc 
comes from observations of quasi-periodic oscillations (QPOs) 
in different accreting stars: neutron stars, black holes, and 



white dwarfs. For example, in low-mass X-ray binaries 
(LMXBs), typical QPO frequencies vary from very high, u ^ 
1300Hz, down to very low, u ^ O.lHz (e.g., [van der Klis| 
[2006 ). An important feature of the spectra is the presence 
of 'twin peaks' of high-frequency QPOs, where the separation 
between peaks is often equal to either the frequency of the 
star or half of its value (e.g., [van der Klis|2000| ). Different the- 
ories have been proposed to explain this phenomenon (e.g., 
Lamb et all|1985[ jMiller, Lamb & Psaltis1|1998[ [Lovelace &| 



Romanova|2007l ). 



In accreting magnetized stars, the magnetosphere may 
rotate slower than the inner disc (e.g., Romanova et al.|2002[ 
[Long et al.[|2005| ). Consequently, there is a maximum in the 
angular velocity distribution as a function of radius at the 
inner disc. This condition can lead to formation of unstable 
radially trapped Rossby waves (Lovela ce & Romanova|200'7l 
[Lovelace eF al. 2009 ). These waves can be responsible for one 
of the high-frequency QPO peaks observed in LMXBs. These 
waves may play a role of blobs suggested earlier in theoreti- 
cal models of QPOs, such as in the beat-frequency model (e.g., 
[Miller, Lamb & Psaltis|1998l [Lamb & Millerf 2003 ) . However, 
compared with the blobs, the inner density waves are much 
more ordered, and have higher coherence. 

Different QPO frequencies in discs around stars and black 
holes can be connected with resonances in the disc, where the 
amplitude of waves may be enhanced (e.g., Kato et al. 1998^ 
Kato.2004i [Zhang & L^2006j. F or example. Lamb & Miller] 
( [2003[ ) and [Lai & Zhang| pOOS) ) suggest that wave can be 
excited at the spin-orbital resonance, and the Hght from the 
high-frequency QPO can interact with this wave and result in 
the lower-frequency QPO. |Kato[ ( [2004| ) suggested that in the 
case of the tilted black holes the disc may be warped, and 
investigated resonance associated with non-linear interaction 
of different waves with this warp (see also |Petri|2005l [Lai &[ 
[Zhang 2008; Kato|2008||2010 ). Slowly-precessing warps can 
form in both the black hole and the magnetic star cases and 
can be responsible for low-frequency QPOs (e.g., |Lai[|1999| 
[Pfeiffer &Lai|2004j ). 

The very low-frequency oscillations (^ 0.1 — IHz) in 
LMXBs can be connected with disc oscillations at large dis- 
tances form the star (e.g., corrugation waves can be excited, 
[Kato et al.||1998| ). |Titarchuk & QsherovichI ( ^200 0) suggested 
that normal (bending g) oscillations of the whole disc may be 
responsible for low-frequency oscillations in accreting com- 
pact stars. Such global oscillations have been observed in nu- 
merical simulations of accretion to a black hole with a tilted 
spin axis (Fragile & Anninos|[2005l [Fragile et al.|[20 07V The 
authors note that the frequency of these oscillations depends 
strongly on the size of the disc. 

In all above-mentioned cases, where a rotating mag- 
netized star excites waves in a disc, the waves have been 
analyzed using linearized equations. The linear theory pro- 
vides an important theoretical understanding of the different 
waves. However, the approach has significant limitations, par- 
ticularly in those cases where the amplitude of waves is not 
small. Further, the theory gives only a rough treatment of the 
force on a disc due to a tilted stellar magnetic field. In this pa- 
per, we present the first results of global three-dimensional 
(3D) MHD simulations of waves excited by a rotating star 
with a misaligned dipole magnetic field. 

In Sec. |2]we briefly describe the linear theory of waves 
in discs, our numerical model and methods for wave analysis 
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in simulations. We also discuss the theory of waves in Keple- 
rian discs in Appendix |a| and the theory of trapped waves in 
Appendix [b] In Sec. [s] we present the results of our simula- 
tions. We discuss applications of the numerical results in Sec. 
|4] Conclusions are given in Sec.js] 



2 THE MODEL, AND ANALYSIS OF WAVES IN THE DISC 

Below we give the highlights of the linear theory of waves 
(Sec. |2.1| ), describe our numerical model (Sec. |2.2| ) and the 
approach which has been used for analysis of waves in nu- 
merical simulations (Sec.|2.3|). 



2 . 1 Theoretical background 

Here we briefly summarize the theory of waves excited in 
accretion discs (see, e.g., Kato et aL]p!998| ). In the majority 
of theoretical studies, the disc is approximated to be thin, 
isothermal and non-magnetic. Furthermore, oscillations can 
be free, or can be excited by the external force. Below we 
briefly discuss different types of waves. 



2.1.1 Free oscillations in the disc 

Small perturbations of a thin, non-magnetized disc involve 
in-plane and out-of-plane displacements of the disc matter. 
Waves can be roughly divided into inertial-acoustic waves 
(referred to as p-modes, where the restoring forces are rota- 
tion and pressure gradient), and inertial- gravity waves (re- 
ferred to as g-modes where the restoring force is gravity). 
Inertial-acoustic modes are in-plane modes, and they prop- 
agate in the form of pressure or density disturbances. They 
are often termed "density waves" (see also Goldreich|1978|). 
The inertial- gravity modes are out-of-plane modes ( Kato et al^ 
[1998 ). These waves lead to bending of the disc, and hence we 
call them "bending waves". 

For both the in-plane and out-of-plane modes, the az- 
imuthal dependence is exp(im(/)), with m — 0, 1, 2, ... The 
m — I in-plane perturbation is a one-armed spiral wave and 
the m — 2 in-plane mode is a two-armed spiral wave. The 
out-of-plane modes involve a small vertical displacement or 
shift of the midplane of the disc, Zw(t, r, 0) (cylindrical coor- 
dinates) with \zw \ <^ r. 

To investigate the propagation of small-amplitude waves 
in the disc, perturbations of the hydro dynamical values Q — 
p, yo, Vr,.- are expanded as (e.g., Okazaki et al.|1987[|T^naka| 
let al.|2002t|Lai & Zhang|2008l ) : 



where ip = 6p/po is the perturbation of the enthalpy for den- 
sity waves or Zu, for the bending waves. Also in this equation, 
u) = (jj — mQ is the Doppler-shifted angular frequency of the 
wave seen by an observer orbiting at the disc's angular ve- 
locity Q(r), Q± is the frequency of oscillations perpendicular 
to the disc, = r~^d{r^Q.^)/dr is the square of the radial 
epicyclic frequency, and Cs is the sound speed in the disc. 

The local free wave solution in the WKBJ approximation 
can be considered as (Kato et al.|1998| ). 



(X k'^ 



-1/2 



exp 



/r 
drkr{r)'\ 



(3) 



for {krvY > 1 and {krhY <^ 1, where kr is the ra- 
dial wavenumber. Thus, equation ([2| can be presented as 
d'^ip/dr'^ = —krip and one gets the dispersion relation: 



k^Cg . (4) 

Propagation of the perturbations is possible in regions where 
(cj^ - k,'^){Cj'^ - nQ\) > 0. The WKBJ approximation breaks 
down in the vicinity of points where the coefficient near ip in 
equation ^ changes sign or has a singularity. At the points 
where uj = ±Hi there is an outer (+sign) or an inner (-sign) 
Lindblad resonance. At the points where uj = ±nQ± there 
are vertical resonances. For cj = there is a corotation reso- 
nance. Appendix|A|summarizes the behavior of these different 
modes. 

In this work we are mainly interested in (1) oscillations 
in the plane of the disc (n — 0), and (2) out of plane os- 
cillations (n = 1). Further, we consider only the one-armed 
(m = 1) and two-armed (m = 2) modes. If the accretion disc 
is Keplerian, then n = = Q. From the analysis of |Zhang| 
|& Lai| f2006), one-armed density waves (m = 1) at n = 
can propagate only in the region of the disc outside the outer 
Lindblad resonance. The inner Lindblad and corotational res- 
onances in this case (n = 0) are absent. For vertical oscilla- 
tions (n = 1), the position of the outer Lindblad and vertical 
resonances coincide. Two-armed density waves (m = 2) with 
n = can propagate either in the external parts of the disc 
outside the outer Lindblad resonance or in the region inside 
the inner Lindblad resonance. 



2.1.2 Waves excited by the external force. Resonances 

In our model, waves are 'driven' by the rotating tilted mag- 
netosphere of the star. The interaction of the magnetosphere 
with the disc is a complex, non-linear, and non-stationary pro- 
cess. For simplicity of the analysis, this interaction can be ap- 
proximated by the force acting on the disk and presented as a 
sum of different harmonics: 



Q(r, z, 0, t) 



^ Qm,nAr)Hn exp(zm</. - tcot) . (1) / = ^' <^ - = E /'"(^' ^) e^P[^ 



(5) 



Here, cj is the angular frequency of the wave, m = 0, 1, 2... 
is the number of arms in the azimuthal direction, Hn,n > 
are the Hermite polynomials (with n = 0, 1, 2, ..), and h{r) 
is the half- thickness of the disc at radius r. Waves with n = 
and n = 1 represent density and bending waves, respectively. 

Linearization of the equations of motion (e.g., |Kato et al.| 
|1998| ) leads to 



where is the angular frequency of the star. In linear ap- 
proximation, we suggest that the m— harmonic of this force 
excites the m— armed wave in the disc. The frequency of this 
force is mQ^ . Then, for a Keplerian disc, the condition for the 
Lindblad resonances uj = ±Q will be mQ* — mQ = or 



TLR 



1 ± 



1 



2/3 



(6) 



dr'^ 



(2) 



where Vcr — {GM is the corotation radius. For bend- 
ing waves, n — 1, the condition for vertical resonances 
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m = 1 m = 2 





^CR 


^OLR 


^OVR 


Hlr 


^OLR 


nvR 


^OVR 






n = 


n = l 


n = 


n = 


n = l 


n = l 


0.544 


1.5 


2.38 


2.38 


0.94 


1.97 


0.94 


1.97 


0.414 


1.8 


2.86 


2.86 


1.13 


2.36 


1.13 


2.36 


0.192 


3.0 


4.76 


4.76 


1.89 


3.93 


1.89 


3.93 


0.089 


5.0 


7.93 


7.93 


3.15 


6.55 


3.15 


6.55 



Table 1. Radii of resonances in a thin hydrodynamic Keplerian disc excited by a rotating perturbation from a magnetized star rotating with 
angular velocity (in dimensionless units discussed in the text). rcR is the radius of the corotation resonance. The next two columns show the 
radii of the outer Lindblad roLR and outer Vertical rovR resonances in the case of one-armed (m = 1) density (n = 0) and bending (n = 1) 
waves. The right four columns show the radii for resonances in the cases of two-arm waves: inner and outer Lindblad resonances tjlr and roLR 
for density waves (n = 0), and inner and outer vertical resonances rjvR and rovR for bending waves (n = 1). 
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/i 


''"cor 







a 






-Rout / -R* 


Warp Properties 


FW/iO.5 


0.5 


1.5 


0.54 


30° 


a 




0.02 


34.5 


fastj 


FW/il.5 


1.5 


1.8 


0.41 


30° 


a 




0.02 


34.5 


fast, Qw ^ 


FW6>5 


0.5 


1.5 


0.54 


5° 


a 




0.02 


34.5 


fast, 


FW6>15 


0.5 


1.5 


0.54 


15° 


a 




0.02 


34.5 


fast, Qw ^ 


FW6>45 


0.5 


1.5 


0.54 


45° 


a 




0.02 


34.5 


fast, 


FW6>60 


0.5 


1.5 


0.54 


60° 


a 




0.02 


34.5 


fast, Qw ^ 


FW6>90 


0.5 


1.5 


0.54 


90° 


a 




0.02 


34.5 


fast, i~ 


FWaO.OO 


0.5 


1.5 


0.54 


30° 


a 




0.00 


34.5 


fast, Qw ^ 


FWa0.04 


0.5 


1.5 


0.54 


30° 


a 




0.04 


34.5 


fast, i~ ^^>i< 


FWa0.06 


0.5 


1.5 


0.54 


30° 


a 




0.06 


34.5 


fast, Qw ^ 


FWa0.08 


0.5 


1.5 


0.54 


30° 


a 




0.08 


34.5 


fast, i~ 


SWcorl.8 


0.5 


1.8 


0.41 


30° 


a 




0.02 


34.5 


slow, Qw ~ 0.517* 


SWcor3 


0.5 


3.0 


0.19 


30° 


a 




0.02 


34.5 


slow, ~ O.IQ^ 


SWcor5 


0.5 


5.0 


0.09 


30° 


a 




0.02 


34.5 


slow, Qw ~ 0.0217* 


LRcorl.8 


0.5 


1.8 


0.41 


30° 


a 




0.02 


57.7 


no big warp 


LRcor3 


0.5 


3.0 


0.19 


30° 


a 




0.02 


57.7 


no big warp 


LRcor5 


0.5 


5.0 


0.09 


30° 


a 




0.02 


57.7 


no big warp 


SRcor3 


0.5 


3.0 


0.19 


30° 


a 




0.02 


26.9 


disc oscillations, example 



Table 2. The main set of simulation runs. See Sec. |3.1| for a detailed description. 



uj = ±Q is the same, and hence the vertical resonances are 
located at the same radii. 

For a one-armed density wave (m = 1, n = 0) in a Kep- 
lerian disc, the outer Lindblad resonance is located at radius 
roLR = {4:GM/Qiy^^. The location of the outer vertical res- 
onance rovR for a bending wave (n = 1) is the same. For a 
two-armed spiral wave (m = 2), the outer Lindblad and verti- 
cal resonances are located at the distance of tolr = ^ovr = 
[9G'M/(4Q^)]^/^ while the inner Lindblad resonance is at the 
distance hlr = nvR = [GM/(4Ql)]^^^, The corotation reso- 
nance for all waves is located at radius rcR = {GM/Qiy^^. 
A more detailed description of the modes is given in Appendix 

m 

In this paper we investigate waves excited by stars with 
different angular velocities, . Table [T] shows the positions 
of corotation, vertical and Lindblad resonances for different 
values of 17*. 



2.1.3 High-frequency trapped waves near magnetized stars 

In relativistic discs, the epicyclic frequency k has a maximum 
at the inner disc (where relativity effects from a black hole 
or a neutron star are stronger), and this may lead to forma- 
tion of waves which are trapped below this maximum (e.g., 
"Okazaki et al. 1987J |Nowak & Wago ner|[T99H [Kluzniak &| 
Abramowicz 2002 ). These waves are important because they 
may be responsible for the high-frequency QPOs observed in 
black hole-hosting systems (e.g., |Stella & Vietri|1999t|Alpar&| 
[Psalti s 2005 , 2008 ). 

In cases where a star has a dynamically-important mag- 
netic field and is slowly rotating, the magnetosphere slows 
down the rotation of matter at the inner part of the disc. 
Hence, there is a maximum in the disc angular velocity dis- 
tribution. This situation is similar to that of relativistic discs: 
high-frequency waves can be trapped below the maximum of 
the angular velocity distribution (e.g., [Lovelace & Romanova] 
|2007[ [Lovelace et al. 12009) ) . We briefly summarize the theory 
of trapped waves around magnetized stars in Appendix |b] 
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Figure 2. Views of the inner parts of the disc in a model with a relatively large magnetic moment, Jl = 1.5 (model FW/il.5) at different times, t. 
Top row: A three-dimensional view of the disc at one density level, p = 0.07. Middle row: z— averaged height of the disc center of mass above and 
below the equatorial plane, Zw Bottom row: surface density distribution in the equatorial plane, a. Time t is given in units of Pq /4, where Pq is 
the orbital period at r = 1. Plots are shown in a coordinate system rotating with the star. 
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Figure 3. Views of bending and density waves shown in Fig.[2]at time t = 80. The left two panels show 3D views of the warp from the side and 
face-on orientations. The right two panels show bending and density waves and the locations of resonances: rcR is the corotation resonance; 
rovR is the outer vertical resonance; tjlr and roLR are the inner and outer Lindblad resonances respectively. 



2.1.4 Low -frequency bending waves 



rotating magnetosphere of the star, in particular if the star 
rotates slowly. 



In the absence of an external force, perturbation of the disc 
in the vertical direction leads to the formation of bending 
waves. The restoring force is the gravity force which is pro- 
portional to —gh/r. The frequency of matter oscillation in 
the vertical direction is Q± ^ Q = y^GM/r^, therefore, at 
large distances from the star, these oscillations have a very 
low frequency. A one-arm global bending wave with low fre- 
quency and long wavelength in the radial direction can prop- 
agate to very large radial distances ( Kato et al. 199 8 ). In addi- 
tion, global bending oscillations of the whole disc are possible 
dTitarchuk & Osherovich||2000| ) . Waves can also be reflected 
or generated at the outer boundary of the disc ( [Zhang & Lai| 
|2006| ) . Such low- frequency oscillations can be excited by the 



2.1.5 More realistic discs 

Real accretion discs have a finite thickness, and hence the thin 
disc approximation may not be applicable. In addition, the 
disc may have a complex distribution of density, temperature 
and magnetic field. All these factors may determine the for- 
mation and propagation of waves (e.g., |Kato et aL]|1998| [Fu| 
|& Lai|2012| ). 

Waves in discs of finite thickness were investigated in 
few global 3D simulations of accretion onto rotating black 
holes with the tilted spin. Formation of warp and tilted discs 
have been observed in these simulations (e.g., [Nelson & Pa-| 
|paloizou[[2000[ [Fragile & Anninos[ [2005] ) . Different disco- 
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Figure 4. Power spectral density (PSD) obtained for waves in model FW/il.5 inside r < 7 of the simulation region (the entire region is Rout = 
12.1). Left panels: PSD for one-armed, m = 1, (top) and two-armed, m = 2, (bottom) bending waves. Right panels: same, but for density waves. 
The solid red lines show the angular velocity distribution in the disc. Dashed lines show angular frequency of the star and positions of resonances. 



seismic modes (such as trapped g-mode), were observed in 
axisymmetric simulations of accretion on to a non-rotating 
black hole from non-magnetized accretion disc (e.g., |0'Neil| 
|et al.|[2009l ). However, in discs where accretion is driven by 
the magneto-rotational instability (MRI, e.g., Bal bus & Haw-| 
|ley||1991) , no interesting waves, such as trapped relativistic 
g-modes were observed (e.g., [Reynolds & Miller|2009| ) . These 
simulations show that the magnetic field may possibly damp 
waves in the disc. 

In our simulations, the disc has a finite thickness, h/r ^ 
0.1, which approximates a thin disc. We observed in simula- 
tions that the angular velocity is almost Keplerian in major 
part of the disc (excluding the inner parts near the magne- 
tosphere) and hence the approximation for Keplerian disc is 
valid. Some magnetic flux of the magnetosphere is trapped 
inside the disc and hence may influence the positions of res- 
onances. However, this is the region where the angular ve- 
locity is strongly non-Keplerian, and where we can not apply 
the standard theory of waves. Overall, results obtained in our 
model can be compared with results obtained with the linear 
theory of waves in a thin disc, excluding the innermost region 
of the disc, where the theory of trapped waves ( |Lovelace"&| 
|Romanova|2007[ [Lovelace et al.|2009l ) is more appropriate. 



2.2 Numerical model and reference values 

We solve the 3D MHD equations with a Godunov-type code 
in a reference frame rotating with the star, using the "cubed 
sphere" grid. The model has been described earlier in a series 
of previous papers (e.g., [Koldoba et al.|2002[ [Romanova et"aL[ 
[2003[|2004l ). Hence, we will describe it only briefly here. 

Initial conditions. A rotating magnetic star is surrounded 
by an accretion disc and a corona. The disc is cold and 
dense, while the corona is hot and rarefied, and at the ref- 
erence point (the inner edge of the disc in the disc plane), 
Tc = lOOTd, and pc = 0.0 Ipd, where subscripts 'd' and 'c' de- 
note the disc and the corona, which are initially in rotational 



hydrodynamic equilibrium (see e.g., [Romanova et al.[2002[ for 
details) . The disc is relatively thin, with the half-thickness to 
radius ratio h/r ^ 0.1. 

Boundary conditions. At both the inner and outer bound- 
aries, most of the variables Aj are taken to have free bound- 
ary conditions, dAj/dr — 0. On the star (the inner boundary) 
the magnetic field is frozen onto the surface of the star. That 
is, the normal component of the field, Bn, is independent of 
time. The other components of the magnetic field vary. At the 
outer boundary, matter is not permitted to flow into the re- 
gion. The simulation region is usually large enough that the 
disc has enough mass to sustain accretion flow for the dura- 
tion of the simulation. The free boundary conditions on the 
hydrodynamic variables at the stellar surface means that ac- 
creting gas can cross the surface of the star without creating a 
disturbance in the star or the flow. This neglects the complex 
physics of interaction between the accreting gas and the star, 
which is expected to produce a strongly non-stationary shock 
wave in the stellar atmosphere (e.g., [Koldoba et al.|2008| ). 

A "cubed sphere" grid is used in the simulations. The grid 
consists of Nr concentric spheres, where each sphere repre- 
sents an inflated cube. Each of the six sides of the inflated 
cube has an x curvilinear Cartesian grid. The whole grid 
consists of 6 X A^r X A^^ cells. The typical grid used in our 
simulations has A^^ = 61^ angular cells in each block. We 
use a different number of grid cells in the radial direction: 
Nr = 130, 140, 160. 

Reference values. The simulations are performed in di- 
mensionless variables A. To obtain the physical dimensional 
values A, the dimensionless values A should be multiplied by 
the corresponding reference values as A = AAq. These 
reference values include: mass Mo = M*, where is the 
mass of the star; distance Rq — i?*/0.35, where R^ is the 
radius of the star; and velocity vq — {GM / RqY^'^ . The time 
scale is the period of rotation at r = i?o : i^o = 2t^Rq /vq . Ref- 
erence angular velocity is uq — vq/ Rq, reference frequency is 
z/Q — cjo/27r. To better resolve the temporal behavior of the 
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Figure 5. Same as in Fig.|2] but for a smaller magnetosphere, /i = 0.5 (model FW/iO.5). 
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Figure 6. Same as in Fig.js] but for a smaller magnetosphere, /i = 0.5 (model FW/iO.5). 



different waves, we record the simulation data every quarter 
of rotation. In the figures we show time t in units of Po/4. 
The reference magnetic moment is /xo = BqRI, where Bq is 
a reference magnetic field. We introduce the dimensionless 
coefficient Jl = /i//io, which acts to control the size of the 
magnetosphere. Here, ji — B^Rl is magnetic moment of the 
star, B^ is the magnetic field of the star at the magnetic equa- 
tor. From the above, we obtain the reference magnetic field 
Bq = {B^rI)/{rI'J1). The reference density is po = Bl/vl 
and the reference mass accretion rate is Mo = pqvqRq. 

Numerical method and code: In our Godunov scheme, we 
use the solver described by Powell et al. (1999). We use a co- 
ordinate system rotating with the star. We split the magnetic 
field into a vacuum dipole componen t and a comp onent cal- 
culated in equations: B = Bd + ( Tanaka||l99 4). Instead 
of using the full energy equation, we use an equation for the 
entropy, which is reasonable for our problems because we do 
not expect the formation of shocks. We use the ideal equation 
of state for the gas, p — Kp^ . Viscosity is added to the code, 
with a— prescription for the viscosity coefficient ( [Shakura"&| 



|Sunyaev||1973| ). Viscosity supports slow, quasi-steady accre- 
tion of matter towards the star during long time intervals. The 
viscosity is nonzero only inside the disc, that is, above a cer- 
tain threshold density (which is p = 0.2 in our dimensionless 
units). We use small value a — 0.02 in most of the simulation 
runs, but take even smaller, as well as larger values in the test 
runs. 



A typical simulation run lasts about 50Po. This time is 
sufficient to observe and analyze most frequencies in the disc. 
To analyze global, low-frequency modes in the disc, we use a 
smaller simulation region. 



We consider only non-relativistic discs. For discs around 
neutron stars the relativistic effects can be important and can 
significantly modify the waves (e.g., |Nowak & Wagoner|199"H 
|Kato et aL]|199 8 ). However, test simulations of waves per- 
formed in the weakly relativistic case show results similar to 
those of the non-relativistic case. 



8 M. M. Romanova et al 




Figure 7. Power spectral density obtained for waves in model FW/iO.5 inside r < 7 of the simulation region. Left panel: PSD for one-armed, 
m = 1, bending wave. Right panel: PSD for two-armed, m = 2, density waves. 
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29 s 


2.2 ms 



Table 3. Sample values of the physical parameters of different types 
of stars. See Sec. |2.2| for a description. 

2.3 Analysis of waves in the disc 

Here, we describe the methods used to analyze the different 
waves observed in our 3D MHD simulations. 



2.3. 1 Visualization of density and bending waves 

To reveal the position and amplitude of density and bending 
waves, we separate the disc from the corona using one of the 
density levels, po = 0.1 y/pjp^ — 0.01, and then integrate 
through the disc in the z— direction. 

To analyze the density waves, we calculate the surface 
density at each point of the disc for different moments of time: 

cr = cr(t, r, 0) = J pdz , (7) 

where r is the radius in the equatorial plane and is the 
azimuthal angle. As a result, we obtain the two-dimensional 
distribution of density waves. 

To analyze the bending waves, we calculate the vertical 
displacement of the center of mass for each point of the disc: 

z^(t,r,0) = ^^-^ , (8) 

and obtain a two-dimensional distribution of bending waves 
for different times. 

We chose a very small value for po in order to be sure 
that we do not miss the regions between spiral waves, or other 
locations with low density. Typically most of matter in the disk 
is concentrated inside the boundary p = 0.1. 

We use the time-dependent distribution of cr(t, r, 0) and 
Zw (t, r, 0) for visuaHzation of density and bending waves, and 
also for spectral analysis described in the next section. 



2.3.2 Spectral analysis of disc oscillations 

Here, we analyze the distribution of the angular frequencies 
of waves in the disc. We consider the variables Zw (t, r, 0) and 
(j(t,r,(j)) and denote them as fj{t,r,(j)), j = 1, 2. This is a 
function of time t, and the two spatial coordinates, r (radius 
in the equatorial plane) and (azimuthal angle). Our simu- 
lations use a coordinate system rotating with the star, and the 
functions / are calculated in this reference frame. At a fixed 
position (r, 0) in the equatorial plane the function / is quasi- 
periodic, because there are always some inhomogeneities in 
the disc which move faster or slower than the coordinate sys- 
tem. They cycle around the rotation axis of the disc. The func- 
tion / is determined at some interval of time, ti < t < t2. It is 
not periodic, but rather oscillates around some value, which 
by itself also varies in time. In addition, it rotates relative to 
a distant observer with the angular velocity of the star, Q*. 
Our objective is to identify the most significant frequencies as 
seen by a distant observer. Also, we want to determine the 
distribution of waves in the disc. 

As a first step of analysis, we "clean" the function /. The 
goal of this cleaning is to exclude from consideration parts of 
the data which are not needed for frequency analysis, for ex- 
ample, the average density, or regular variations of the func- 
tion /. In general /(ti) / f{t2), and in the observer's frame, 
/ r^'(t2) . We consider the function which has been 
obtained as a result of this cleaning as a function determined 
for the whole time-interval t2 — ti, with period T — t2 — ti. 
In order to do this, we subtract from /(t) the linear function: 

12 — ti 

and then subtract from the resulting function the time- 
averaged value 

fit) ^ fit) - (/). 

In the next step, we select the azimuthal modes of the func- 
tion /. The function /°^^ (which is in the observer's frame) is 
connected to function / in the star's frame by the relationship: 

r^^(t,r, 0) = /(t,r,(/)-Q*t). 

We use the Fourier expansion of /°^^: 

/°^^(t, r, 0) = ^ /m ^(t, r) exp(im(/)) , 

where 

= — / d(f)f{t,r,(f) — QU)e^p{—im(f)) 

27T J 





= (ix/(t,r,x)exp[-im(x+r2*t)] = /,n(t,r)exp(-im^7*t). 

In the next step, we analyze the function = 
fm exp(— imQ*t). Functions fm and are also determined 
in the interval ti <t<t2- Subsequently, we perform a spec- 
tral analysis of /m, and interpret the term exp(— imll^t) as 
the displacement of the frequency by the value mQ*. 

Next, we perform the time Fourier transform of /m(t, r). 
We consider fm as a periodic function of time with period 
T — t2 — ti. In this case, the angular frequency has a dis- 
crete set of values with step 5uj — 27r/T. The coefficients of 
the Fourier transformation are (here we neglect the numbers 
corresponding to frequencies) : 

= i y dtfm{t,r) e:>cp{iujt-imQU) - fm{oo-mQ.^) - fm{oj)- 

Here, Co = uj — mQ*. Hence, the Fourier coefficients of the 
function fm, determined in the rotating coordinate system, 
are calculated for frequencies considered in the rotating coor- 
dinate system. When we transfer to the observer's coordinate 
system, all frequencies are misplaced: cj = cD + mQ*, but the 
Fourier coefficients do not change. 

To characterize the spectral properties of the function 
/m(t,r), we use the value of the Power Spectral Density 
(PSD): 

PSD=\f^%uj,r)\\ 
This value characterizes the power in the signal with fre- 



quency cj. Here, we analyze one-armed (m = 1) and two- 
armed (m = 2) modes for density and bending waves. 



3 WAVES OBSERVED IN 3D MHD SIMULATIONS 

We performed multiple simulation runs for different parame- 
ters of the star at fixed parameters of the disc (see Sec. |3.1| 
and Table [2j. The results strongly depend on the relative 
positions of the magnetospheric radius, rm We can sepa- 
rate our results into two main groups: (1) The inner disc ro- 
tates with angular velocity close to that of the star (that is, 
Vm ~ Vcr)' In this case, the rotating magnetosphere excites 
a high-amplitude bending wave (a warp), which rotates with 
the angular velocity of the star (see Sec. |3.2| ). (2) The inner 
disc rotates faster than the star (that is, rm < Vcr)- In these 
models, the high-amplitude (corotating with the star) warp 
does not form (see Sec. |3.3| ). 

3.1 Description of simulation runs 

We use as a base case a model with the following parameters: 
the magnetic moment Ji — 0.5 (in dimensionless units, see 
Sec. |2.2| ), so that the magnetospheric radius obtained from 
simulations is rm ~ 1.3 (or, in stellar radii, Vm ~ 3.7i^*). 
Such a magnetosphere is sufficiently large to excite waves in 
the disc. At the same time the simulation runs have a rea- 
sonable duration (on a supercomputer running at 168 pro- 
cessors), so that we were able to run multiple cases with 

^ rm is a radius, where the matter stress in the disc matches the 
magnetic stress in the magnetosphere and the corotation radius, rcr- 
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Figure 9. Analysis of bending waves in the disc in cases where the magnetic moment of the star is tilted at angles ^ = 5°, 15° and 30° (models 
FW^5, FW^15, and FW^30). Top panels: xz-slices of density distribution and sample magnetic field lines. Middle panels: 3D view of matter flow 
at one density level, p = 0.07. Bottom panels: the height Zw of the disc center mass. Dashed circles show the location of corotation and vertical 
resonances: rcR = 1.5, rovR = 2.38. 




different parameters. The star rotates w^ith angular velocity 
= 0.54 (in dimensionless units). The corresponding coro- 
tation radius is rcr = 1-5 is only slightly (1.15 times) larger 
than the magnetospheric radius. We take the angle betw^een 
the star's magnetic moment and the spin axis (the misalign- 
ment angle) to be ^ = 30°. The maximum (outer) radius of 
the simulation region is Rout = 12.1 (in our dimensionless 
units), or Rout = 12.1/0.35 = 34. 5i?* (in radii of the star). 
We use the described model as a base and call it FW/xO.S. Ab- 
breviation "FW" indicates a "Fast Warp" because in all cases 
where rm ~ Vcr, a strong bending wave (a warp) forms and 



rotates with the angular velocity of the star, which is fast, 
compared to other cases. 

We varied different parameters, taking model FW/xG.S as 
a base model. In model FW/il.5 we show an example of a 
larger magnetosphere, Jl — 1.5, where similar waves are ex- 
cited but with larger amplitudes (see Sec. 3.2.1|). We discuss 



the warp amplitude in this model in Sec. |3.2.3 



In models FW6'5, FW6'15, FW6'45, FW6>60 and FW6'90 we 
investigate excitation of waves for different misalignment an- 
gles and compare them with the base case FW/iO.5, which 
corresponds to 6* = 30° (see Sec. |3.2.4| ). 

In models FWaO.OO, FWa0.04, FWa0.06, FWa0.06 and 
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Figure 11. Warps observed in simulations where the viscosity parameter a varies from a = 0.0 to o; = 0.08 (models FWaO.OO - FWaO.OS ). 



FWaO.08, we investigate the possible dependence of warp 
properties on the a— coefficient of viscosity (see Sec. |3.2.5| . 

In another set of simulation runs we drop the condition 
Vm ~ Vcr and decrease the angular velocity of the star, 
(thereby increasing the corotation radius rcr). We observed 
that the strong magnetospheric warp disappears but other 
t3^es of waves become stronger. In some cases, a slow warp 
is observed (see models SWcorl.8, SWcorS and SWcorS). Ex- 
periments with different sizes of the simulation region show 
that this warp appears only in models with relatively small 
simulation region, where free disc oscillations are excited 
more rapidly. We demonstrate excitation of free disc oscilla- 
tions using a smaller simulation region. Rout — 9.4 = 26. 9i?* 
(model SRcorS). In another set of simulation runs with a 
large region. Rout — 20.2 = 57. 7 R^, a warp is not present 
(models LRcorl.8, LRcorS and LRcorS). However, apart from 
the main magnetospheric warp, which can be present or not, 
other types of bending and density waves were observed in 
all simulation runs. 



3.2 Waves and Fast Warps in Cases where rm ~ Vcr 

Here, we investigate the case of rapidly-rotating stars, where 
the corotation radius is approximately equal to the magneto- 
spheric radius, rm ~ Vcr- This case is interesting because it 
approximately corresponds to the rotational equilibrium state 
which is the most probable state in the life of accreting mag- 
netized stars (e.g., [Long et al.|2005t ). 

We performed two simulation runs (models FW/xO.S and 
FW/il.5) for stars with different sizes of the magnetosphere 
(fi = 0.5 and Jl — 1.5) and matching corotation radii rcr — 
1.5 and rcr = 1.8 (so that rm ~ ^cr). The dipole magneto- 
sphere is tilted at ^ = 30°. We observed that in both cases, a 
large warp forms at the corotation radius, and it rotates with 
the angular velocity of the star. We use these two cases to 
demonstrate the formation of magnetospheric warps and to 
show their main properties. 



3.2. 1 Formation of fast warp and other waves in the case of a 
large magnetosphere (Model FWiil.5) 

First, we investigate the case of large magnetospheres (model 
FW/il.5), where the warp is larger than in the base case. We 
observed from simulations that the magnetosphere truncated 
the disc at the distance of rm ~ 1.55 from the star. This ra- 
dius rm is 1.16 times smaller than rcr — 1.8 and therefore 
rm ~ Tcr- Fig. [T] shows that matter of the inner disc has been 
lifted above the disc and flows to the star in two ordered fun- 
nel streams (see [Romanova et al.|2003[[2004| ). It can also be 



seen that the inner parts of the disc have a bent structure. In 
particular, the far left regions of the disc are lifted above the 
rest of the disc, forming a warp. 

Fig. [2] (top panels) shows three-dimensional views of the 
inner disc density distribution (one of the density levels) at 
different moments in time (plots are shown in the coordi- 
nate system rotating with the star). One can see that the 
high-amplitude bending wave (a warp), seen in Fig.[T] is lo- 
cated at approximately the same position in time, and thus 
it corotates with the star and its magnetosphere. The middle 
panels show the height of the center of mass of the disc Zu, 
(see eq. [sj above and below the equatorial plane. It can be 
seen that the location of the warp coincides with the inner, 
high-amplitude part of the bending wave, which also coro- 
tates with the star (red color). The bending wave is evident 
on the opposite side of the disc (blue color). Thus, we can 
see that the bending wave represents a one-armed {m — 1) 
feature which propagates to larger distances. The bottom pan- 
els show the averaged density distribution in the accretion 
disc, cr (see eq. [t]) . One can see that the two-armed (m = 2) 
trailing density wave (n = 0) forms at the corotation radius 
and propagates to larger distances. The formation of the warp 
which corotates with the magnetosphere and breaks at some 
distance from the star is in agreement with theoretical predic- 
tions ( [Terquem & Papalo izou 2000|. The warp does not pre- 
cess (as discussed by |Lainl999j ), because in our simulations 
the rotational axis of the dipole is aligned with the rotational 
axis of the disc (like in [Terquem & Papaloizou|2000l ) . 

The left two panels in Fig. [s] show an expanded view of 
the warp shown in Fig. [2] both in side and face projections. 
The right two panels show the locations of resonances calcu- 
lated theoretically for a thin Keplerian disc (see Table [T] and 
Appendix|A]). One can see that the bending wave starts at the 
corotation radius. It passes through the Outer Vertical Res- 
onance (OVR) and propagates farther out. According to the 
theory (e.g., K ato et al.|1998|rzhang & Lai >2006, see Sec.jAl] 
of the Appendix A) a bending wave is expected to have its 
maximum amplitude at and around vertical resonance (see 
Fig.|A2]of the Appendix, 3rd panel from the left). Fig.[3]shows 
that in our simulations, a warp really has the maximum am- 
plitude around OVR, which is only slightly shifted towards 
CR. The shape of the bending wave observed in our simula- 
tions is similar to that derived from the theoretical analysis 
(see right panel of Fig.|A2]in Sec. |Al| ). The right-hand panel 
shows the distribution of density waves, and the locations of 
the resonances. 

According to the theory (e.g., [Zhang & Lai|[2006| ), one- 
armed (m = 1) in-plane (n = 0) density waves can propagate 
to the exterior of OLR, while two-armed (m = 2) waves can 
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Figure 12. Waves observed in models with large simulation region {Rout = 20.2 = 57. 7 R^) at different corotation radii. From top to bottom: 
models LRcorl.8, LRcorS, and LRcorS. Left panels: 3D view at one of the density levels (p); Middle panels: bending waves (zt^— distribution); 
Right panels: density waves (cj— distribution). 



propagate to the exterior of OLR and to the interior of ILR. 
In our case, there are clear, two-armed density w^aves vN^hich 
propagate at r > tolr, which is in agreement with the theory. 
The ILR is located inside the magnetosphere and therefore, 
there is no disc at r < tilr. It is possible that these density 
waves are excited directly by the magnetic forcing on the in- 
ner parts of the disc. Alternatively, the magnetic force may 
excite predominantly bending waves, which in turn compress 
matter and generate density waves. A one-to-one comparison 
between bending and density waves shows some similarity in 
the pattern between them, so that this is a possibility. 

We apply the spectral analysis described in Sec. |2.3.2| to 
this model (FW/il.5). The top left panel of Fig. [4| shows PSD 
for the one-armed (m = 1) bending (n = 1) waves. It can 
be seen that PSD has a clear feature corresponding to a warp, 
which rotates rigidly with the angular velocity of the star. This 



feature starts at the corotation radius and extends outward to 
a distance ofr ^ 2.7rm- There are also lower-frequency bend- 
ing waves, which are at frequencies above and below Keple- 
rian. The higher-frequency waves may represent a coupling 
between Keplerian vertical oscillations and bending waves ex- 
cited by the star. Waves of the lowest frequency may represent 
a coupling between Keplerian frequency and the frequency 
of free, bending oscillations of the disc. Analysis of the two- 
armed (m = 2) bending waves shows that these waves are 
much weaker with a PSD negligibly small compared with the 
one-armed bending wave. 

We also calculated the PSD for one-armed and two- 
armed density waves, and observed that the PSD of a two- 
armed wave is much stronger. Fig. [4| (bottom right panel) 
shows that for radii tolr < r < 5, there is a feature with 
angular frequency u ^ 0.25 — 0.3 which corresponds to the 
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density waves observed in Fig. [2] (bottom panels). This fre- 
quency corresponds to the frequency of both waves. In addi- 
tion, PSD is high for several frequencies at r < rcR. One of 
these frequencies (the highest) may correspond to trapped 
density waves (e.g., [Lovelace et al.| [2009] ) (see Sec. |3.3.2| ). 
These waves appear where the angular velocity of the disc 
decreases to match the angular velocity of the star for the 
case where the magnetosphere rotates slower than the inner 
disc. In this model, the corotation radius is only slightly larger 
than the magnetospheric radius. Nevertheless, there is a max- 
imum in the angular velocity distribution (solid red line in 
Fig. [4]) where trapped waves can form. Fig. |4] (right bottom 
panel) shows a candidate trapped wave, which is located in- 
terior to the peak in angular velocity distribution at r ^ 1.5 
and its angular velocity a; = 0.9 which is slightly higher than 
2Q* ^0.83. There is another density wave at about the same 
distance, r ^ 1.5 — 1.6 but at a lower frequency, ^ 0.7 — 0.8. 
This may be a wave excited by magnetic force at the disc- 
magnetosphere boundary, but it weakens at the corotation ra- 
dius. There is also a wave at frequency 2Q^ with smaller PSD. 
This wave may be excited by the warp in the disc. One-armed 
density waves have negligibly small PSD (right top panel) . 

It is interesting to note that the PSD for density waves 
is very different from that for bending waves. This supports 
the first hypothesis discussed above that bending and density 
waves are excited independently by the tilted magnetosphere. 

3.2.2 Formation of fast warp and other waves in the case of a 
smaller magnetosphere (FWfiO.S - base model) 

We calculated another similar case, but for a star with a 
smaller magnetic moment, jl — 0.5. In this case the disc is 
truncated at a smaller radius, rm ~ 1-2. We took a smaller 
corotation radius, Vcr — 1-5, in order to have Vm ~ rcr 
as in model FW/il.5 (here, we obtain sHghtly larger value, 
Tcr/rm ~ 1.25 compared with the model FW/il.5, where 
Tcr/rm ~ 1.15). We observed that this case is similar to model 
FW/il.5, however all scales are smaller because of the smaller 
size of the magnetosphere, which excites the waves. 

Fig. [5] (top panels) show a 3D view of bending waves 
and the inner warp. The middle panels show that a bending 
wave forms at the corotation radius and propagates outward. 
The inner part of the bending wave (a warp) is located at the 
same position indicating that the warp corotates with the star 
The bottom panels show the formation of two-armed density 
waves, similar to model FW//1.5, but at smaller scale. 

Fig. [6] (left two panels) shows expanded views of the 
warp in two projections. The two right panels show the po- 
sitions of bending and density waves relative to resonances. 
One can see that the warp part of the bending wave is located 
between the corotation and the vertical resonances, with max- 
imum amplitude at the vertical resonance. The right-hand 
panel shows that density waves are excited at the Lindblad 
resonance but are damped rapidly at r ^ 3.4, which approx- 
imately corresponds to the radial extent of the warp. In ad- 
dition, one can see density enhancements between the closed 
magnetosphere and the corotation radius. These are trapped 
density waves. 

Fig. [7] (left panel) shows that in the PSD plot the main 
feature corresponds to the warp which rotates with the angu- 
lar velocity of the star and spreads out to r 3, which is 2.5 
radii of closed magnetosphere. 



The right-hand panel of Fig. [t] shows that the two-armed 
outer density wave is located at tolr < r < 4. The angu- 
lar frequency (for both waves) is uj ^ 0.4. PSD also shows 
a two-armed trapped density wave with angular frequency 
uj ^ 1.14 which is located at the inner disc, with r ^ 1.15. 
The frequency of this wave is located below the maximum in 
the disc angular velocity distribution. 

Thus, the warp and other features in model FW/iO.5 are 
similar to those of model FW/il.5, which has a larger mag- 
netosphere. In both cases, the warp forms near the magneto- 
sphere and corotates with the star The warp extends out to 
r ^ (2.5 — 2.7) rm and rotates rigidly in this region. Below, we 
use a model with a larger magnetosphere, FW/il.5, to deter- 
mine the amplitude of the warp. In all other models, we use 
model FW/iO.5 as a base case. For this model, the simulation 
runs are not as long, and we can investigate the model at a 
number of different parameters (see Table [2j. 

3.2.3 Amplitude of the warp 

It is important to know the ampHtude of the warp. A large- 
amplitude warp can block or diminish the light from a star 
In fact, this mechanism of light-blocking by a warped inner 
disc has been proposed as an explanations for dips in the 
Hght-curves observed in some T Tauri stars ( [Bouvier et al.| 
|2003l |2007b[ [Alencar et al.|[2010| . Clearly the warp ampli- 
tude needs to be sufficiently large to produce the observed 
dips. 

Figures [2] and [s] (middle panels) show the amplitude of 
the warp, Zw (the position of the center of mass of the disc rel- 
ative to the equatorial plane) . Warp amplitude is the largest, 
reaching Zw ^ 0.2, in model FW/il.5. Taking into account the 
distance of this maximum to the star, r ^ 2,we obtain the ra- 
tio: Zw/r ^ 0.1. In the model with a smaller magnetosphere, 
FW/iO.5, we find ^ 0.1 - 0.15, and the ratio z^/r ^ 0.05. 
These numbers are in agreement with the theoretical anal- 
ysis of warped discs by Terquem & P apaloizou| ( |2000| ), who 
concluded that the amplitude of the warp is expected to be 
5 — 10%. However, they considered the case of a very thin 
disc. 

In a realistic situation, the disc has a finite thickness, and 
the height of the warp hw that can obscure light from the star 
can be larger than the deviation of the center of mass from 
the equatorial plane, Zw The density in the disc decreases 
with height, and the height of the warp also depends on the 
chosen density level, hw = hwip)- 

To determine the height of the warp from our simula- 
tions, we chose model FW/il.5, which has a larger magneto- 
sphere. Fig. [s] (bottom panel) shows a view of the magneto- 
sphere and the inner warp, while the top panel shows a slice 
of the density distribution in the xz plane, which crosses the 
warp (note that the slice does not go through the maximum 
amplitude of the warp, which is at about 30° — 40° clockwise 
from the xz— plane). 

The top panel shows that the height of the warp depends 
on the chosen density level: the lower the density level, the 
higher the warp. In the figures, we chose a set of density lev- 
els. One can see that at the highest density levels (red color), 
the warp amplitude is not very large, while at the lowest 
shown density level (dark blue) the warp rises up to 30% of 
the distance to the star ihw/r = 0.3) and can be significant in 
blocking the light at different inclination angles of observer 
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The warp amplitude is even larger at even lower density lev- 
els. However, for an analysis of light obscuration by the warp, 
it is important to know the opacity at different density levels 
in each particular case. We note that the stellar light can also 
be obscured by the funnel streams, which are located closer 
to the star than the warp, and can provide an even larger ef- 
ficient height of the obscuring material, hw/r > 0.3. 

The relative amplitude z^/hu, and the height of the 
warp, hw/r, depend on the force which magnetosphere ap- 
plies to the inner disc. It is expected that both values should 
increase with the size of the magnetosphere, that is with fl. 
Comparison of amplitudes done in previous sections, shows 
that the ratio Zw/r is indeed twice as large in case of larger 
magnetosphere, jl = 1.5, than in small magnetosphere, (jl = 
0.5). 

It may also depend on the temperature of the inner disc 
and its thickness. Fig. [s] shows that the accretion disc is thin, 
compared with the height of the warp and hence the magnetic 
force is the main one which determines the height of the warp 
in this model. 

We observed from our simulations that the ratios z^/r 
and hw/r are larger in model FW/xl.S (with a larger magneto- 
sphere, rm ~ 4.7i?*), compared with model FW/iO.5 (where 
the magnetospheric radius is smaller, rm ~ 3.3i?*). This is be- 
cause a larger magnetosphere appHes stronger magnetic force 
to the disc (which has the same properties in these two mod- 
els) . We suggest, that an even stronger warp can be expected 
in the cases of an even larger magnetosphere. It is expected 
that in many CTTSs, the disc can be truncated at larger dis- 
tances (e.g., Tm ~ (7 - 8)i?* in AA Tau star, Donati 2011| ). 
Special investigations are needed to understand whether the 
amplitude and height of the warp will increase or decrease 
with a further increase of rm . 



3.2.4 Warps in the cases of different tilts of the dipole, 

Here, we investigate the dependence of the warp amplitude 
Zw on the misalignment angle of the dipole moment relative 
to the star's spin axis, 6. For this reason, we performed a series 
of simulations at different 6 values, from ^ = 5° up to ^ = 
90° (models FW6'5, FW6'15, FW6'45, FW6'60, and FW6'90). We 
take as a base model FW/iO.5, where 6 — 30°. 

Fig. [9] shows the results for cases of relatively small tilts 
of the dipole: = 5°, 15°, and 30°. The top panels show 
sHces of density distribution in the xz— plane and sample field 
lines. One can see that the inner disc is thicker in models with 
smaller 6. The densest part of the disc (red color) shows bend- 
ing waves. The middle panels show a 3D view of the bending 
waves at one density level. One can see that waves are excited 
in all three models. The bottom panels show the amplitude of 
the warp, Zw . A warp forms in all three cases, with the max- 
imum at the vertical resonance (outer dashed circle). From 
these plots, we find the values Zw/r ^ 0.025, 0.045, and 0.05 
in models with 6* = 5°, 15°, and 30°, respectively The height 
of the warp, hw is larger than Zw . Top panels of Fig. [9] show 
that at small values of ^ (5° and 15°), matter accumulates at 
the disc-magnetosphere boundary. The disc thickens, and the 
inner disc can absorb stellar Hght, thus providing hw/r ^ 0.3. 
M. 6 — 30°, the amplitude of the warp Zw is larger than in 
the other two cases. However, the height of the inner disc is 
smaller. 

The bottom panels of Fig. [9] show that in cases where 



6' = 5° and 15°, some matter accretes to the star through a 
Raleigh-Taylor instability (e.g., [Romanova et al.|2008t|Kulka-| 
|rni & R omanova 2008 ) . However, this does not prevent excita- 
tion of bending waves, because these phenomena occur at dif- 
ferent radii: unstable accretion occurs at the inner regions of 
the disc, while bending waves are excited at somewhat larger 
distances by the external parts of the magnetosphere. 

Fig. [To] shows the excitation of bending waves around 
stars with larger misalignment angles: 6 — 45°, 60° and 90°. 
Top two panels show that bending waves are excited in all 
of these cases. However, the bottom panel shows that a high- 
amplitude warp forms only when — 45° and 60° (with ratios 
z^/r — 0.04 and Zw/r — 0.03, respectively), while in the case 
of ^ = 90°, where the magnetosphere is almost symmetric 
(some of the non-axisymmetry is connected with the initial 
distribution of magnetic flux in the disc), ratio Zw/r is about 
10 — 20 times smaller than in the other two cases. 



3.2.5 Influence of viscosity on the warps 

In our code an a— type viscosity is incorporated into the disc. 
Viscosity is thus regulated by parameter a (see [Shakura "&| 
|Sunyaev|1973| . All of the above-mentioned simulations were 
performed for a relatively small viscosity, a — 0.02. Here, 
we show the results of our simulations for the base model 
FW/iO.5 but with different viscosities: a = 0,0.04,0.06 and 
0.08. 

Fig. [TT] shows bending waves for different a values. One 
can see that in all of the cases a warp forms, and is located 
at approximately the same place, with a maximum amplitude 
between corotation and vertical resonances, or at the vertical 
resonance. Here, for larger values of a, we show the plot for 
a later moment in time, t, because in this case the accretion 
rate is larger, the disc has higher density, and bending waves 
are excited more slowly, compared with the case of a smaller 
a. 

The possible influence of viscosity on the evolution of 
warped discs has been investigated earlier (Bardee n & Pet- 
19751 iPapaloizou & Pringle||1983t Papaloizou & Lin 



terson 



1995| lTerquem & Papaloizou^2000V6gilvie 2006). Papaloizou 



& Pringle ( 1983 ) (see also Papaloizou & Lin 1995 ) considered 



the formation of warps in viscous a— type discs, which have a 
vertical scale height h(r), and concluded that to a first approx- 
imation, the warp satisfies a wave-type equation if a < h/r. 
In the opposite case, if a > h/r, the warp satisfies a diffusion- 
type equation. 

In our simulations, the ratio h/r ^ 0.1. Hence, we ex- 
pect non-diffusive, wave-like behavior for a values smaller 
than a — 0.1. This may explain why in all the cases consid- 
ered above, where a < 0.08, we observed a warp with similar 
properties. 



3.3 Waves in cases where the magnetosphere rotates 
slower than the inner disc, rm < rcr 

In the above section we investigated cases where the mag- 
netosphere and the inner disc rotate approximately with the 
same angular velocity, which corresponds to the rotational 
equilibrium state (e.g., [Long et al.|[200 5!). We observed that 
the main feature is a warp corotating with the star. How- 
ever, the situation is different when the magnetosphere ro- 
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Figure 15. Same as in Fig.[T3] but for model LRcorS (large region, rcr = 5). 
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tates much slower than the inner disc, that is, when the coro- 
tation radius rcr is much larger than the magnetospheric ra- 
dius, Tm. Such a situation may arise during periods of en- 
hanced mass accretion, when the disc compresses the mag- 
netosphere and the magnetospheric radius rm decreases. Al- 
ternatively, the dipole component of the magnetic field of a 
young star may decrease due to internal dynamo processes 
(e.g., |Donati 2011). As a result, the magnetospheric radius 
rm will decrease, providing the condition rm < rcr- 

If the magnetosphere rotates slower than the inner disc, 
then it applies force to the inner disc, but with relatively low 
frequency of the star's rotation, Q*. This will excite bending 
and density waves with the frequency of the star. On the other 
hand, the matter of the inner disc has a much higher angular 
velocity. It interacts with the slowly-rotating magnetosphere, 
which serves as a non-axisymmetric obstacle for the disc mat- 
ter flow. It exerts a non-axisymmetric force on the inner parts 
of the disc, thus exciting higher-frequency waves (both den- 
sity and bending waves) . Thus, a complex pattern of waves is 
expected in this case. In addition, we observed that a slowly- 
rotating star excites global bending oscillations of the whole 
disc if the size of the simulation region (that is, the size of the 
disc) is not very large. In the case of a large simulation region, 
these waves are also excited, but their frequency is very low 
and they are excited slowly compared with the overall time 
of simulations. This is why we initially consider the case of a 
relatively large simulation region, where these low-frequency 
waves can be neglected (see Sec. |3.3.1| ). Next, we investi- 
gate the cases of smaller simulation regions, where global disc 
oscillations develop and couple with other frequency modes 
(see Sec.[3l4l). 

3.3.1 Investigation of waves in a large disc 

Here, we investigate the waves excited in a relatively large 
simulation region, with radius Rmax = 20.2 ^ 57. 7i?*. This 
region is 1.7 times larger than that used in Sec. |3.2[ We ob- 
serve that the very low-frequency global mode of disc oscilla- 
tions develops slowly, and we are able to investigate higher- 
frequency oscillations at the "background" of this slowly- 
growing mode. 

We use a model with parameters corresponding to our 
base model FW/iO.5, but take a larger simulation region and a 
larger corotation radius: rcr = 1.8, 3 and 5 (models LRcorl.8, 
LRcorS, and LRcorS). The magnetospheric radius rm ~ 1-15 
is located at the same place as in model FW/xO.S. Hence, the 
new models correspond to stars with slower rotation. The 
simulations show that matter in the rapidly- rotating inner disc 
interacts with the slowly-rotating tilted magnetosphere of the 
star, so that different types of waves are excited and propa- 
gate to larger distances. 

Fig. [12] (top row) shows the results of our simulations for 
model LRcorl.8, where the corotation radius is rcr — 1.8, and 
the star rotates relatively rapidly compared with the other two 
cases. The left and middle panels show that a bending wave is 
excited at the disc-magnetosphere boundary and propagates 
outward, forming a spiral wave. The amplitude is enhanced 
at the outer vertical resonance, r = rovR (see Table [l]), and 
this enhancement has a position similar to that of the warp in 
the base case (FW/iO.5), where the wave corotates with the 
star. Here, however, instead of corotation, we observe a ten- 
dency to corotate, where a wave rotates with the frequency 



of the star during only a part of the rotational phase, then it 
slows down, but soon after, another similar wave appears and 
corotates with the star, and so on. The right panel shows sur- 
face density distribution, a. One can see that different waves 
form at different radii from the magnetosphere, including a 
two-armed density wave, which forms at the outer Lindblad 
resonance. 

Next, we show the results for model LRcorS, where the 
corotation radius is larger, rcr — 3. Fig. [12] (middle row)) 
shows that bending waves are less ordered inside the corota- 
tion radius compared with model LRcorl.8. They look like 
parts of a tightly wound spiral wave extending out to the 
corotation radius. At larger distances, a new bending wave 
is formed, and it is more ordered and not as tightly wound. 
This observation is in agreement with the theoretical predic- 
tion, where a tightly wound spiral wave is expected at r < rcr 
(see the left two panels of Fig. |A2] in AppendixjA]), and a much 
less tightly wound spiral is expected at r > rcr (see Fig. |A2| 
right panels) . The right-hand panel of Fig. [12] (middle row) 
shows that different density waves formed between the inner 
radius of the disc and the OLR radius. 

Next, we show the results for model LRcorS, where the 
corotation radius is even larger, rcr — 5. Fig. [12] (bottom 
row) shows that bending waves are not ordered and have a 
rippled structure, as in the rcr — 3 case. The middle panel 
shows that bending waves are not ordered, but one can track 
a tightly wound spiral wave inside the corotation radius (in 
agreement with Fig.|A2]of Appendix |a]). A more ordered, but 
less tightly wound spiral wave starts dXr ^ rcr and becomes 
yet more ordered at the OVR. The right-hand panel shows 
density waves. One can see that different density waves prop- 
agate at r < riLR' There are waves at larger distances, but 
they are not seen at the chosen density levels. 

The PSD shows interesting sets of bending and density 
waves in these three models (see Figures [13) [T4] and [T5| . The 
PSD confirms that there is no large-scale warp at the stel- 
lar frequency, which is different from models FW/xO.S and 
FW/il.5 (see Sec. |3.2| ). Instead, only a small-scale bending 
wave is excited at the stellar frequency or at twice this value 
(see left panels of the figures for m = 1 and 2) . This wave is 
only present at the inner region of the disc, close to the region 
where matter lifted to the funnel flow, and may therefore re- 
flect this lifting. For the three models considered, the angular 
frequency of the wave is low. For example, in the case of a 
one-armed wave (m = 1), cj ^ 0.41, 0.19 and 0.09 for models 
LRcorl.8, LRcorS and LRcorS, respectively. 

There are also high-frequency bending waves associated 
with interaction of the rapidly-rotating matter of the inner 
disc with the non-axisymmetric magnetosphere. The angular 
frequency of these waves, cj 0.8 — 1, approximately corre- 
sponds to Keplerian velocity of the inner disc. Also, there are 
bending waves of intermediate frequency. All of these waves 
are located at the inner edge of the disc, 0.8 < r < 1.2 where 
matter has already been lifted, or has started lifting to the 
funnel flow. 

Lower-frequency bending waves excited by the slowly- 
rotating star propagate from the disc-magnetosphere bound- 
ary to large distances, and we often see an enhanced PSD at 
the stellar frequency, and at different distances from the star. 
The PSD can be further enhanced at the corotation radius, as 
we can see in Fig. [l4] (top left panel) . The stellar- frequency 
mode is also seen in other models, but at lower values of PSD. 
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The density waves show a variety of frequencies (see 
right-hand panels of figures [Ts] [l4| and [T5| . At low frequen- 
cies, we observe enhanced PSD at radii ^3 — 4, which often 
correlate with a position of different resonance (e.g., in model 
LRcorl.8 with OLR, in model LRcorS with CR and OLR, in 
model LRcorS - with ILR). Plots of the surface density dis- 
tribution a for these models (see Fig. [T2] ) show that in all 
these models, the density waves have higher amplitudes at 
the distances of r < 4, where the influence of the rotating 
magnetosphere is stronger. 



3.3.2 Waves in the inner disc: trapped density waves 

At high frequencies, the PSD of density waves is often en- 
hanced at radii corresponding to the peak in the angular ve- 
locity distribution, and the frequencies of these waves are usu- 
ally somewhat lower than the peak frequency (see Fig. [T3][T4| 
and [is] where the solid red lines show the angular velocity 
distribution in the disc). We suggest that these are trapped 
density waves, which are expected in cases where the angu- 
lar velocity distribution in the disc has a maximum ( [Lovelace] 
|& Romanova|2007[ [Lovelace et al.|2009D . In our models, the 
maximum in the angular velocity as a function of r occurs 
because the star rotates slower than the inner disc. The re- 
gion where the angular velocity decreases gives rise to the 
possibility of radially trapped unstable Rossby waves with az- 
imuthal mode numbers m > 1 and radial width Ar/r <^ 1 
( [Lovelace et al.[[20Q9 ). T he theory of trapped waves is sum- 
marized in Appendix [Bj Here, we analyze trapped waves in 
one of our models (LRcorS), where these waves are clearly 
observed. Fig. [T4| shows a high-frequency (cj ^ 0.57), one- 
armed (m = 1) wave, which is located below the maximum 
in uj at radii 1 < r < 1.5. To check this wave we also plot 
the surface density distribution at the inner disc for different 
times and clearly see that there is a one-armed wave located 
approximately between radii 1 and 1.5 (see Fig.[l6]). We be- 
lieve that this is a trapped density wave. According to the 
theory ( ^Lovelace et al. 2009 ) , trapped waves have a tendency 
to be within the maximum in the u distribution. We can see 
that in model LRcorS, the maximum of the PSD for m = 1 is 
located in the left part of the cj— curve, while for model LR- 
corS, the trapped wave is located below the maximum of cj 
for both the m = 1 and m — 2 modes. 

The bottom panels of Fig. [16] show bending waves for 
the same moments in time. A small inner warp forms at 
1 < r < 1.5, but rotates with the angular velocity of the star. 
This small warp is also seen in Fig. [l4] (top left panel) . The 
PSD also shows the presence of a high-frequency, (cj ^ 1), 
bending wave at r < 1 (top left panel). This feature prob- 
ably corresponds to a small funnel or tongue, which forms 
due to the magnetic Raleigh-Taylor (RT) instability at the disc- 
magnetosphere boundary. 



3.3.3 Accretion through instabilities: density and bending 
waves 

At the disc-magnetosphere boundary, matter can penetrate 
through the magnetosphere due to the magnetic Rayleigh- 
Taylor (RT) instability (e .g., |Arons & Lea|1976[[Spruit & Taam] 
[1990[ [Kaisig et al.|1992| ) . A necessary condition for the insta- 
bility is that the gravitational acceleration towards the star be 



larger than the centrifugal acceleration. Additional factors are 
also important, such as the gradient of the angular velocity in 
the inner disc and the gradient of the surface density (e.g., 
[Spruit et al.|1995j . Global numerical simulations of accretion 
onto a star with a tilted dipole magnetic field show that accre- 



tion may be in the stable or unstable regimes (see Romanova 
|et al.[[2008[ [Kulkarni & Romanova|[2008[ [205 9: BachettTet 
al. 2010). In the stable regime, matter accretes in ordered 
funnel streams, which corotate with the star. In the unstable 
regime, matter accretes through unstable temporary funnels 
or tongues, which carry the angular momentum of the inner 
disc, and can rotate faster than the inner disc matter. The 
simulations show that the criterion for the onset of the RT 
instability is close to that derived byl Spruit et al.[ ( [T995i ). In 
addition, the boundary between stable and unstable regimes 
depends on the tilt of the magnetosphere and on the simula- 
tion grid[^ At larger tilt angles, more matter accretes in stable 
funnel streams. 

In our models with rcr = 3 and 5, a star rotates slowly 
relative to the inner disc, and therefore the instability is fa- 
vored. On the other hand, the tilt of the dipole is relatively 
large {_0 = 30°), and this favors accretion through funnel 
streams. This is why in the case of rcr — 3, where the star 
rotates relatively slowly, we see only a weak instability. To 
demonstrate accretion through instabilities, we take an even 
slower rotating star, where Vcr — 5 (model LRcorS) and ac- 
cretion through the RT instability is very clear. The top panels 
of Fig. [it] show that the instability leads to the formation of 
one main "tongue", which rotates with the angular frequency 
u ^ 0.5. This unstable tongue is associated with the main 
one-armed trapped density wave that forms in the inner disc. 
The frequency of this trapped wave is clearly seen in Fig. [Ts] 
(top right panel). The PSD also shows the presence of a sec- 
ond harmonic (m — 2) with cj ^ 1.1 at radius r ^ 1 (see 
bottom right panel of Fig. [T5|. This wave is barely seen in the 
surface density plots in Fig.[17[ but is clearly seen in the PSD. 

The bottom panels of Fig. [Tt] show that the instability is 
also seen in the bending waves that form at the inner edge of 
the disc. The two-armed bending wave (bottom left panel of 
Fig. [Ts]) also has a mode cj ^ 1.15 at r < 1, which is close 
to the mode cj = 1.1 of the m = 2 trapped density wave. 
Here, we have an example of how a trapped density wave 
generates unstable tongues at smaller radii. Another possible 
explanation is that both waves indicate accretion through in- 
stabilities. 

Fig. [is] (top left panel) shows the presence of a high- 
frequency (cj ^ 0.9 — 1.1), one-armed (m = 1) bending wave 
at r < 1. A similar mode is observed in two other models (LR- 
corl.8, LRcorS). This mode is probably connected with the 
rotation of unstable funnel streams at the inner edge of the 
disc. The frequency of the mode is higher than that of the disc 
(solid red line) . We suggest that this mode may reflect the ro- 
tation of temporary unstable funnel streams, which are lifted 
above the magnetosphere, and therefore do not experience 
magnetic braking (as matter in the disc does), and hence can 



^ Recent simulations performed with a finer grid (as in this paper) 
show that the RT instability occurs more readily compared with the 



earlier simulations performed with a coarser grid (e.g., , Romanova et[ 
[al.|2008| . 
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rotate with a nearly Keplerian velocity. We call these waves 
inner bending waves. 

3.4 Waves in a smaller simulation region: global disc 
oscillations and slowly-rotating warps 

In sections [3. 3.1 13.3.3) we considered the case of a relatively 
large region, where free bending oscillations in the disc de- 
velop slowly. Here, we took a smaller simulation region and 
observed that these oscillations develop more rapidly (see, 
e.g., |Titarchuk & Osherovich,200Q| for a possible mechanism 
of global bending oscillations of the disc) . 

First, we demonstrate the excitation of free bending oscil- 
lations in the disc of a very small size. Rout — 9.4 = 26. 9i?* 
(model SRcor3). The simulations show that initially, waves 
are excited at the disc-magnetosphere boundary (see top pan- 
els of Fig.[l8]). Later, at t 70 - 100, global bending oscilla- 
tions of the disc develop. However, at t > 100, the maximum 
of the bending wave amplitude moves to the region between 
the OVR and the CR resonances (see Fig. [Ts] ) . This feature is 
similar to the warps in cases of rapidly-rotating stars (see Sec. 
|3.2| ) . However, this new warp rotates slower than the star. 

Next, we investigate the waves in the intermediate-sized 
simulation region. Rout — 12.1 = 34. 5i?*, where free bend- 
ing oscillations in the disc develop sufficiently fast in particu- 
lar in cases of slowly-rotating stars) . We consider two models 
with corotation radii Vcr = 1.8 and 3 (models SWcorl.8 and 
SWcor3). 



In model SWcorl.8, we observed evolution of bending 
waves similar to that in model SRcor3 (see Fig. [Ts] ) . As a re- 
sult of this evolution, a warp forms at radii rcR < r < rovR. 
This warp rotates slower than the star (see Fig. [T9] ) . A closer 
view of the warped disc (see Fig. [2Q| ) shows that a signifi- 
cant part of the inner disc is tilted and involved in this slow 
rotation. 

In the case of an even more slowly-rotating star, rcr — 3 
(model SWcor3), a very slow warp forms at r > rcR. Fig. [21] 
shows the rotation of this warp about the star. 

Fig. [22] shows the PSD for these two models. One can see 
that in the case where Vcr — 1-8 case (left panel), the warp 
rotates with an angular frequency u ^ 0.1, which is about 
four times smaller than frequency of the star. The PSD also 
shows that the main part of the warp is located between the 
corotation and the vertical resonances. In the case of a slower 
rotating star (rcr = 3, right panel), the warp rotates very 
slowly, with an angular frequency cj = 0.01 — 0.02, though 
the accuracy of our model is not very high in cases of very 
low frequencies (because the number of rotations per simula- 
tion run is small) . One can see that in both cases, the angular 
frequency of slow warps is much smaller than the local Kep- 
lerian frequency. In model SWcor3, the angular frequency of 
the warp is comparable with the Keplerian frequency at the 
edge of the simulation region, (joxir ^ 12) ^ 0.024. How- 
ever, this is probably just a coincidence. In test simulations 
with even more slowly- rotating stars (rcr = 5), a warp forms, 
and its angular frequency is even lower: cj < 0.01. We did not 
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establish how this slow warp forms. Its frequency probably 
represents a coupling between the frequency of the star and 
the frequency of global oscillations in the disc. 



4 APPLICATIONS AND DISCUSSION 

The results of our simulations can be applied to various types 
of stars where the magnetospheric radius is several times 
larger than the radius of the star|^ including, for example, 
classical T Tauri stars, accreting millisecond pulsars that host 
spun-up neutron stars and some types of cataclysmic variables 
(such as Dwarf Novae) . Possible applications are briefly dis- 
cussed in sections below. 

4.1 Obscuration of light by warped discs in classical T 
Tauri stars (CTTS) and drifting period 

The light curves of classical T Tauri stars (CTTS) display 
photometric, spectroscopic and polarimetric variations on 
timescales ranging from a few hours to several weeks ( Herbst 
|eFal..2QQ2) . |Bouvier et al.| ( [19991 ) studied the variability of 
AA Tau in great detail and concluded that the photometric 
variabihty with a period of 8.2 days which is comparable to 
the expected rotation period of the star is due to the occul- 
tation of the star by a warp formed in the inner disc (the 
system is observed almost edge-on) . They proposed that this 
warp is produced by the interaction of the disc with the stellar 
magnetic dipole tilted with respect to the disc rotational axis. 
Later, more detailed analysis confirmed the hypothesis of ob- 
scuration by the disc ( [Bouvier et al.||2003[ |2QQ7b). Doppler 
tomography observations of AA Tau show that the dominant 
component of the field is a 2 — 3 kG dipole field that is tilted 
at 20° relative to the rotational axis. At this field strength 
the magnetospheric radius is close to the corotation radius: 
Vm ~ Vcr dPonati et al. 12010) ) . Our simulations show that in 
this situation a large amplitude warp forms and rotates with 

^ Modeling of much larger magnetospheres with truncation radii 
Tm » lOR*, requires much longer simulation runs and is at present 
not practical. Special efforts are required to investigate stars with 
large magnetospheres, such as X-ray pulsars. 



the frequency of the star, and the AA Tau warp may be similar 
to that described in model FW/il.5. An analysis of the warps 
in cases where the dipole has different tilts (see Sec. |3.2.4| ) 
shows that an angle of ^ < 20° is sufficiently large to generate 
a warp. It is possible that warps are common features in many 
CTTSs. For example, [Alencar et aL] ( [2010| ) analyzed the pho- 
tometric variability of CTTSs in the young cluster NGC 2264 
using data obtained by the CoRoT satellite and concluded that 
AA Tau-like light curves are fairly common, and are present 
in at least 30 — 40% of young stars with inner dusty discs (see 
also [Carpenter et al.|2001| ) . 

Not only fast warps, but also other waves can be gener- 
ated by the rotating magnetosphere in CTTSs. However, it is 
difficult to observe these waves in CTTSs: these stars are usu- 
ally brighter than the inner disc and are strongly variable due 
to coronal magnetic activity. In addition, observations usually 
cover not more than a few periods of the star, and it is difficult 
to extract a possible frequency of the waves from light-curves. 
One possibility is that the waves generated in the inner disc 
may influence the formation of magnetospheric streams and 
their subsequent motion, thereby leaving a trace on the sur- 
face of the star in the form of hot spots, which can move faster 
or slower than the star ( [Romanova et al.|2004l [Bachetti et al.| 
|2010| ). Observations show that many CTTSs do not have a 
precise period, but instead a quasi-period, which varies with 
time around some value (e.g., [Rucinski et al.||2008| ). This 
quasi-period may be connected with the formation of high- 
frequency waves, such as trapped density waves. These waves 
represent regions of enhanced density and therefore probable 
places where funnel streams may form. The position and fre- 
quency of a trapped waves vary with accretion rate. Thus, the 
angular velocity of hot spots on the surface of the star and the 
associated quasi-period of the star will also vary. 

In a number of situations accretion through the RT 
instability is possible, where matter accretes in unstable, 
stochastically-forming tongues (e.g., [Romanova et al.|[2"008t 
'Kulkarni & Romano va|2008| ) . In this case, the light-curve be- 
comes less ordered, and in a strongly unstable case a CTTS 
may show the absence of a period. It is often the case that 
the instability is not very strong, and both stochastic and pe- 
riodic components are expected in the frequency spectrum 
( [Kulkarni & Romanova|2009T ) . In these cases, trapped density 
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Figure 19. A slowly-rotating warp observed in model SWcorl.8. Top row: 3D views of the warp at different moments in time; Bottom row: center 
of mass of the disc Zw, which demonstrates a slowly-rotating warp. 




Figure 20. An example of a warp which rotates slower than the star in model SWcorl.8 (simulation region Rout = 12.1 = 34. 5i?*). One of the 
density levels is shown. Vector Q^w shows the direction perpendicular to the inner disc warp. 



waves can determine the position of unstable tongues. In the 
case of small magnetospheres, this leads to the formation of 
two coherent tongues that rotate more rapidly than the star 
( [Romanova & Ku lkarni'2009) . Hence, one of the frequencies 
can be connected with trapped density waves. On the other 
hand, we often observe in the PSD even higher- frequency os- 
cillations (inner bending waves) that have a nearly Keplerian 
frequency. Both of these waves may leave an imprint at the 
surface of CTTSs in the form of moving hot spots. Our simu- 
lations show that the coherence of trapped waves is probably 
higher than that of inner bending waves, and it is more likely 
that the variable quasi-period of CTTSs is connected with the 
trapped waves. 



4.2 Application to Accreting Millisecond Pulsars 

Accreting millisecond pulsars (AMPs) show various types of 
quasi-periodic oscillations (QPOs), ranging from very high 
frequencies (^ 1,300 Hz) down to very low frequencies 
(^ O.lHz) (e.g. van der Klis| |2006). One of the prominent 
features in the spectrum of AMPs is a pair of QPO peaks with 
an upper frequency Vu and a lower frequency vi, which move 
in pairs. The distance between the peaks Uu — i^i usually cor- 
responds to either the frequency of the star, z/*, to half this 
value, zy*/2 (e.g., |van der Klis|200()l|2006|), or, in some AMPs 



there is no such correlation (e.g., Mendez e t al.|1998 



Bout- 



loukos et a"Ll|2006[ [Belloni et aiyl ]2007; M endez & Belloni 
2007| ). There are also QPOs with lower frequencies, including 
hecto-hertz QPOs with a frequency of Uh ~ i^u/^, as well as 
low frequency QPOs with a frequency of z/lf ~ z^^^/lS. There 



are also very low- frequency oscillations, vll < IHz. Below 
we discuss possible connections between the waves observed 
in simulations and some of these frequencies. 

4.2. 1 High-frequency oscillations in the inner disc 

Our simulations show that there are two main types of high- 
frequency waves: (1) trapped density waves, and (2) inner 
bending waves. The frequency of the inner bending waves 
is almost Keplerian and is higher than the frequency of the 
trapped density waves, which is located below the peak in the 
disc angular velocity distribution. We suggest that the trapped 
waves can be responsible for the lower QPO peak, vi, and the 
inner bending wave for the upper QPO peak, Uu - 

The frequencies of both waves are associated with the 
inner disc and therefore both will vary with accretion rate 
and will move in pairs. We expect that the quality factor 
Q — vj of both waves will increase with accretion rate, 
because when the disc moves closer to the star, the maxi- 
mum in the angular velocity distribution becomes more pro- 
nounced, and the trapped waves become stronger and more 
coherent. At the same time, accretion through instabilities be- 
comes stronger and more ordered (where one or two tongues 
become dominant and create a more ordered motion, see 
[Romanova & Kulkarni||2009j ) . This is in agreement with the 
observations of QPOs, where Q increases with frequency for 
both upper and lower QPOs (e.g., [van der Klis 2000 ; Barret] 
|et al.|2006|| ) . The quality factor of lower QPOs starts decreas- 
ing at very high frequencies, possibly due to relativistic effects 
(e.g., |Barret et al.|2006l ). 
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Waves in the inner disc can influence the magnetospheric 
flow and can therefore leave imprints on the frequencies of 
moving spots on the surface of the star. Moving spots on the 
surface of the star were observed in global 3D simulations 
of stable magnetospheric accretion ([Romanova~et al. 2003 , 



2004|) and during accre tion through instabilities (e.g., |Kulka-| 
rni & Romanova|2008l ). |Bachetti et al.| ( |2010l ) found that ro- 



tating spots provide two QPO frequencies: one corresponds 
to rotating funnel streams and has a lower frequency, while 
the other corresponds to rotating unstable tongues and has a 
higher frequency (see also [Kulkarni & Romanova||2009| ) . We 
suggest that these two frequencies represent an imprint of the 
trapped and inner bending waves, respectively 

We can calculate the difference in frequency between 
these two high-frequency waves in the disc. For example, in 
model LRcor3, the frequency of the m — 1 trapped wave is 
cJtrap ~ 0.55, while the frequency of the inner bending wave 
is oohend ~ 0.95 (see top panels in Fig. [14]). Then, using Table 
[3] for millisecond pulsars, we obtain the dimensional frequen- 
cies of these waves, vi — 247Hz and Vu — 427¥Lz, and the 
difference between them, Az/ = 180Hz. Similarly, in model 
LRcorS (see top panels of Fig. [Ts]) we obtain utrap ~ 0.51 
(ui ^ 229Hz) for the trapped density waves, and two candi- 
date frequencies for the inner bending waves, Ubend ~ 0.91 
and iObend ~ 0.99, which correspond to i^u = 409Hz and 
jy^ = 445Hz. The difference between the frequencies is 180Hz 
and 21 GHz, respectively. This difference is typical for many 
accreting millisecond pulsars (see bottom panel of fig.l of 
iMendez & Belloni|2007l ). 

For case of moving spots, we investigated how frequen- 



cies vary with accretion rate and found that the typical differ- 
ence is Aiy ^ 200 — 250Hz, which is in agreement with the 
frequency difference for high-frequency waves. This strength- 
ens our hypothesis that the hot spots represent an imprint of 
the waves in the disc. 

We should note that the frequency difference Au is ap- 
proximately the same for stars with rcr = 3 as for stars with 
rcr = 5, and therefore Au does not depend on the frequency 
of the star. This model can explain those AMPs where Au does 
not correlate with the frequency of the star (e.g., Mendez et 



al.|1998[[Boutloukos et al.|2006||B'elloni et al.|2007[|Mend'ez 
& Belloni|2007F ^ 

It has been suggested earlier that the high-frequency 
QPOs may be connected with the rotation of some blobs in 
the inner parts of the disc (Lamb et al.^l985j . The position of 
the inner disc can be determined, for example, by the magne- 
tospheric radius (e.g., [Roman o va & Kulkarni'2009'"Bachet ti et| 
al. 2010 ), or by a sonic point (Miller, Lamb & Psaltis 1998). In 
our simulations, instead of blobs we observe trapped density 
waves or inner bending waves. These waves are more ordered 
than blobs and can provide a much higher quality factor. 



4.2.2 Frequency Modulation at the Vertical Resonance 

In many AMPs, the difference between the QPO peaks is equal 
to the frequency of the star zy* or to z/*/2 (e.g., van der Klis 
|2000| ). Here, we discuss a possible model that follows from 
our simulations and that satisfies to one or both conditions. 

A theoretical analysis of the disc waves shows that if 
waves are excited by the vertical force oscillating with the 
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frequency Vf, then the largest ampHtude of the waves is ex- 
pected at the outer vertical resonance, where the frequency 
of the wave is z^ovr = z^// 2. [Lamb & Miller ( 2003 ) proposed 
a model where the upper frequency, Uu, is connected to the 
rotation of blobs in the inner disc, while the lower frequency 
vi is connected to the modulation of light at this resonance, 
ui = — i/f /2. In the case of a tilted rotating dipole magne- 
tosphere, the forcing frequency can be either one or two times 
the frequency of the star, and therefore the difference in fre 
quencies Uu — i^i is equal to either z/* or z/*/2 (Lai & Zhang 
|2008j ) . This model presents a promising possibility. However, 
the problem is that we do not see a bending wave with the 
predicted frequency at the OVR. 

In our simulations, a warp is usually located at the ver- 
tical resonance (or between OVR and CR) . However, the fre- 
quency of the warp either equals the frequency of the star (if 
rm ~ Tcr, see Sec. |3.2| ), has a low frequency (if the simula- 
tion region is small, see Sec. |3.4| ), or there is no large warp (if 
Tm < Tcr and the simulation region is large, see Sec. |3. 3.1] ). 
From these possibilities we choose the most probable one, 
where rm ~ rcr and the high-amplitude warp corotates with 
the star. 

We suggest that the upper QPO is connected with one 
of the high-frequency waves described in Sec. |4.2.1| (or with 
the rotation of hot spots associated with this wave) . The light 
from this high-frequency QPO can be reprocessed by a warp, 
corotating with the star, so that the expected frequency of the 
lower QPO is ui = z/^i — zy* . Such a model can explain the cases 
where the frequency difference equals the frequency of the 
star. This approach is reminiscent of the beat-frequency model 
(e.g.. Lamb 1995), but in our model, we outline a possible 
mechanism for the "beat" frequency. 

4.2.3 Low -frequency, free bending oscillations in the disc 

Our simulations show that the tilted magnetosphere excites 
different types of waves in the disc, including bending oscilla- 
tions of the whole disc. These oscillations have very low fre- 
quency, which can be of the order of the Keplerian frequency 
at the outer radius of the disc, or even lower. 

If the simulation region is relatively small, then these 
low-frequency waves couple with the higher-frequency bend- 
ing waves excited by the rotating magnetosphere, and a 
slowly-rotating warp forms with the maximum of the ampli- 
tude at the OVR, and with the frequency lower than the Kep- 
lerian frequency at the OVR. 

One of the observed low-frequency QPOs has a frequency 
of ulf, about 15 times lower than that of the high-frequency 
QPOs (e.g., I van der K Hs 2006). We can estimate the size of 
the disc in cases where the frequency is equal to Keplerian 
frequency at the edge of the disc. If the upper-frequency QPO 
is generated at distance ru from the star (it can be, for ex- 
ample, the truncation radius of the disc), then the Keplerian 
frequency at distance r is z/(r) = Vuirulr)^^^ - For the fre- 
quency ratio Vul^LF — 15, we obtain the radius correspond- 
ing to vlf, rLF ^ 6.1r^i. Free oscillations of the inner part of 
the disc with this radius are possible if the properties of the 
disc inside this radius differ from the rest of the disc. We sug- 
gest that the magnetic-to-matter ratio may be larger inside 
this region than at larger distances. For example, if accre- 
tion is provided by the magneto-rotational instabiUty (e.g., 
[Balbus & Hawley||1991| ), then the turbulent cells in the in- 



ner disc become strongly stretched in the azimuthal direc- 
tion, and the azimuthal component of the field is generated 
due to rapid azimuthal rotation (e.g., Hawley|2000[ |Armitage| 
|2002[ [Romanova et al. 2012 ). This is only an example, and 
there may be other reasons for different properties in the in- 
ner disc. There are, of course, other possibilities. For example, 
|Lai| ( [1999 )) suggested that the low-frequency oscillations may 
be connected with a warp precessing around a magnetized 
star, which is expected in cases where the rotational axes of 
the star and the disc are misaUgned (see also [Pfeiffer & Lai| 
[2004] ). 

To explain the even lower-frequency oscillations, such 
as z/LL ^ IHz, we need to consider an even larger disc. 
For example, if the high frequency ~ 10^ Hz) oscilla- 
tions are generated at distance ru, then the radius corre- 
sponding to IHz Keplerian frequency is located at distance 
^LL ^ (yul^iA^Y^'^ru ^ 100r^i, which may be comparable to 
the overall size of the disc. 



4.3 Application to Dwarf Nova oscillations 

The Dwarf Novae (DN) type of cataclysmic variables reveal 
two main types of oscillations: the "Dwarf Nova Oscillations" 
(DNOs) and the "Quasi-Periodic Oscillations" (QPO s) (e.g., 
|Patterson|ri9811 1 Warner] [20041 [Pretorius et al.||2006) . DNO 
oscillations are observed exclusively during outbursts. They 
have relatively short periods, ranging from about 7 to 40 s, 
and have a very high quality factor. They display the period- 
luminosity correlation. 

QPOs usually occur in dwarf novae outbursts, but are 
sometimes seen during the quiescence. They have longer peri- 
ods than the DNO, ranging from about 30 to 1000 s. The QPO 
frequency is usually 15 times lower than the DNO frequency. 
QPOs are characterized by broad humps in the power spec- 
trum, indicating their short coherence time of only a few cy- 
cles. [Robinson & Nather| ( |1979| ) suggested that these oscilla- 
tions are produced in the accretion disc, possibly in a "travel- 
Hng spiral pattern" or in "oscillating rings". Alternatively, they 
can be connected with the precession of a warp ( |Lai||1999| 
[Pfeiffer &La"i|2004 ). 

In our model, we suggest that dwarf novae have a small 
magnetosphere, and the DNO oscillations can be explained 
by trapped density waves, which develop in the inner disc. 
They have a high frequency, which varies with the accretion 
rate (luminosity) . These oscillations are expected to be highly 
coherent. There are a few possibilities for explaining QPOs 
(see Sec. |4.2j ). However, it is difficult to explain the ratio 
T-'DNo /t-'qpo ~ 15, which is the same for various accreting 
compact stars and also for black holes (e.g., [Belloni et al.[ 
[20021 [Warner et al. 12003] ). 



5 CONCLUSIONS 

We performed the first global three-dimensional simulations 
of waves in a disc excited by a rotating magnetized star with 
a tilted dipole magnetic field. The main conclusions are the 
following: 

L In cases where the inner disc approximately corotates 
with the star (rm ~ ^cr): 

1. A strong one-armed bending wave (warp) forms near 
the corotation radius and has its maximum amplitude at the 
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radius of the vertical resonance. It corotates with the magne- 
tosphere and propagates out to the distance of r ^ (2 — 2.5) rm 
from the star. 

2. In the warp, the height of the center of mass of the 
disc, Zw, can reach (5 - 10)% of the distance to the star. This 
is in agreement with values found by |Terquem & Papaloizou| 
( |2000| ) in their study for thin discs. However, the height of 
the warp hw (in our discs of finite thickness) can be as high 
as 30%. 

3. Warps are excited by stars with dipole fields tilted at 
different angles, 0° < ^ < 90°. However, the amplitude of 
the warp is larger in cases of intermediate tilts of the dipole, 
15° < 6* < 60°. 

4. Density waves are excited at different radii from the 
star, but their amplitudes decrease beyond the distance of 
r ^ 4, which may be connected with the relatively large thick- 
ness of the disc (/i/r ^ 0.1). Bending waves propagate in the 
whole simulation region, and have large amplitudes at differ- 
ent distances from the star. 

II. In cases where the magnetosphere rotates more slowly 
than the inner disc (rm < Vcr): 

5. We observe that the strong inner warp is absent. 
Only a much smaller warp corotates with the star. Instead, 
a tightly wound rippled bending wave forms at the disc- 
magnetosphere boundary and propagates out to the corota- 
tion radius. Beyond the corotation radius, the wave is more 
ordered and less tightly wrapped, and it propagates to large 
distances. 

6. High-frequency trapped density waves form in the inner 
disc, within the maximum in the angular velocity distribution. 
These waves were predicted theoretically by [Lovelace & Ro-| 
[manova ( 2007 ) and Lovel ace et aLl C2009). 

7. High-frequency waves of another type - inner bending 
waves - originate at the disc-magnetosphere boundary during 
periods of unstable accretion. 

8. Low- frequency global bending oscillations of the 
whole disc are excited in all cases. In the case of a smaller 
simulation region, we observe a slowly-rotating warp with a 
maximum amplitude at the vertical resonance. The low, sub- 
Keplerian frequency results from the coupling between the 
frequency of the star and the global oscillation mode of the 
disc. In the case of a relatively large simulation region, a 
slowly-rotating warp does not form. Instead, global bending 
oscillations are observed. 

III. Some Applications: 

9. The formation of a high-amplitude warp corotating 
with the star is important for understanding obscuration of 
light in young stars (e.g., |Bouvier et al.|2003t ). The formation 
of warps is expected in other types of stars. The details of ob- 
scuration depend on the transparency of matter in the warp 
to the star's radiation. 

10. Trapped density waves and inner bending waves can be 
responsible for high-frequency QPOs in accreting millisecond 
pulsars, CVs, and for drifting "period" in CTTSs. Trapped den- 
sity waves can appear in a wide variety of situations, where 
the magnetosphere rotates slower than the inner disc. 

1 1 . Inner high-frequency waves can determine the posi- 
tion of funnel streams and moving spots on the surface of the 
star. Hence, they can leave an imprint on QPO frequencies as- 
sociated with moving spots discussed e.g. by |Bachetti et al.| 
( |20T0l ). 



,= l,n=l C 




m=l, n=0 . 



no propagation / OLR 



r/rcr 



b 1 







m=2, n=l 


- 




OVR / 






11 / 


- 


-r* 1— 


1 1 1 1 




Figure Al. Radial wavenumber kr as a function of radius r for a 
thin Keplerian disc with h/r = Cs/{rQ) = 0.05, where rcr is the 
corotation radius. Panels (a) and (b) are for the bending waves, while 
(c) and (d) are for the density waves. 



12. Global disc oscillations and slowly- rotating warps can 
be responsible for low- frequency QPOs in different stars. 
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APPENDIX A: WAVES IN A THIN KEPLERIAN DISC 

Here, we briefly summarize the key aspects of the low mode- 
number waves in thin Keplerian discs. 



Al m = 1 Bending Wave 

For a one-armed (m = 1) bending wave (n = 1) of a Keple- 
rian disc = Q and Q± = Q) equation (4) can be written as 



krCs 



~ \ 2 



1 



0, (Al) 



Q and the other parameters are defined in 



where cj = 
Sec. 2. 

We find that the WKBJ equation (3) is a good approxi- 
mation to the full wave solution except for the region where 
\krrcr \ is small or the wavelength is long. Specifically, the ap- 
proximation is reasonable for \krrcr\ > 4.94. For \krrcr\ < 5, 
equation (Al) is solved as a differential equation for ob- 
tained by setting k^ -d'^/dr'^. 
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Figure A2. From left to right: (1) Profiles of the real and imaginary parts of the displacement of the disc (r, (p = 0) for the 'inner warp', 
r/rcr < 1, assuming Cs/{rQ) = 0.05. (2) Vertical displacement of the disc Zw{r, (f)) for the 'inner warp', r/rcr < 1, assuming Cs/{rQ) = 0.05. 
Here, white = negative and red = positive. (3) Real and imaginary parts of the disc displacement z^: for the 'outer warp' from equation Al. (4) 
Vertical displacement of the disc Zw{r, cj)) for the 'outer warp' for 1 < r/rcr < 2.5 assuming Cs/{rQ) = 0.05. 



Panel (a) of Figure Al shows the radial dependence of 
IkrVcrl for a disc with constant fractional half- thickness h/r — 
Cs/(rQ) — 0.05 and u — ^(rcr). The corotation resonance 
where a; = and /c^ ^ oo is at r = rcr- The vertical/Lindblad 
resonance where Co — and kl = is at r = ~ 
1.59rcr. For r /vcr 0, krVcr 2{r /h)^r /rcr- For r /vcr 

1, krTcr ^ {r /h)\{r/rcr)^^'^ - 1|~\ For r/rcr OO, krVcr 

{r/h)^r/rcr. 

We refer to the warp in the range r /rcr < 1 as the 'inner 
warp.' For most of the range r /vcr < 1 (excluding very small 
r/rcr), the amplitude of the inner warp is well approximated 
by equation (A2). Fig. A2 (left panel) shows the real and 
imaginary parts of the warp profile at = 0. Next panel to 
the right shows the top view of ?R:[zw (r, 0)] of the inner warp. 
The warp becomes very tightly wrapped as r/rcr increases 
towards unity. Approximately, ^ — (2r/3/i) \n[rcr/{rcr — r)] 
and \zw \ ^ \/rcr - r as r ^ rcr. 

In the range r/rcr > we refer to warp as the 'outer 
warp'. Fig. A2 (3rd panel) shows the real and imaginary parts 
of at = 0. As r ^ Vcr-lzw I ^ Vr — rcr while for r > rcr 
\zw\ ^ r~^/^. Fig. A2 (right panel) shows the top view of 



A2 m= 2 Bending Wave 

Using the WKBJ approximation to equation (2) for a two- 
armed (m = 2) warp (n = 1) of a Keplerian disc gives 



krTc 



r 1 

h f 



[(f^/^-l)^-l/4][4-(f^/^-l)-^] 



(A2) 



where f = r /rcr with rcr the corotation radius, so that cj — 
mQ(rcr)- Note that /cr = at the inner and outer vertical 
resonances (IVR and OVR) at f = (1 ± 1/2)^/^ ^ 0.630, 1.31, 
and that \kr \ ^ oo at the corotation resonance at f = 1. Panel 
(b) of Fig. Al shows the radial dependence of kr for a sample 
case. 



A3 In-plane modes, n = (Density Waves) 

The WKBJ solution of equation (2) forn = 0, m = 1, 2, and 
a Keplerian disc gives 



where we have set k = Q. For a disc with Cs/{rQ) = h/r = 
const, one finds 



1/2 



(A4) 



\krrcr\ = {r/h){l/r) [m^r'^^ - if - l] 

where f = r/rcr with rcr being the corotation radius, cj — 
mQ{rcr)- For the one-armed mode (m = 1), there are real 
values of kr only for r/rcr > 2^^^ ^ 1.59, which is the outer 
Lindblad resonance. 

For the two-armed mode (m = 2), there are real values 
of kr only for r/rcr < 2~^/^ ^ 0.63 (the inner Lindblad res- 
onance) and for r/rcr > 1.5^/^ ^ 1.311 (the outer Lindblad 
resonance) . Panels (c) and (d) show the the radial dependen- 
cies of kr for the two cases. 



APPENDIX B: RADIALLY TRAPPED MODES 

A radially localized, in-plane instability of discs for m = 
1, 2, .. know as the 'Rossby wave instability' (RWI) may occur 
near regions where there is an extremum (as a function of r) 
of the entropy and/or the potential vorticity z • ( V x v) /H, 
where v is the disc velocity and S is the disc's surface mass 
de nsity ([Lovelace et aLl|1999| |Li et al.|[2000t [Lovelace et 
ar'2Q09^ 'Lovelace & Romanova||200'7y ; see also [Lovelace 
& Hohlfeld (1978 ). Radially trapped modes in discs around 
black holes were analyzed earlier by Nowa k & Wagoner] 
( ^1991) . Of particular interest here is the RWI in the non- 
Keplerian region of a disc ( [Lovelace et al.[|2009[ hereafter 
LTR09) which occurs when a disc encounter the magneto- 
sphere of a slowly-rotating star. An example of such a rotation 
curve is shown in Figure Bl. For given profiles of Vd), p, etc. 
and a given m, the value of the growth rate is found by solving 
for the 'ground state' solution of Schrodinger-like equation for 
the enthalpy perturbation ip, 



(Bl) 



oj^ = Q"^ -\- ki(?s 



(A3) 



dr'^ 

where U{r\uj) is an effective potential (LTR09) shown for a 
sample case in Figure B2. The ground state is found by vary- 
ing both the real and imaginary parts of the mode frequency 
uj so as to give the V^(r) most deeply trapped in the poten- 
tial well. For typical profiles and m — 1, the growth rate 
is Ui 0.l{vcf)/r)R and the real part of the frequency is 
ujr ~ {vcf)/r)R, where the i?— subscript indicates evaluation 



Warps, Bending and Density Waves Excited by Rotating Magnetized Stars: 25 





\ V 
K 


- bX 




/ : 








3 






1400 
1200 
1000 
800 
600 
400 
200 




Figure Bl. Left panel: Midplane radial dependencies of the azimuthal 
disc velocity magnetic field Bz, and density p, where vk is the 
Keplerian velocity The vertical arrow indicating is the location 
of the bottom of a potential well U{r) shown in right panel (from 
LTR09). Right panel: Effective potential U{r) for the profiles shown 
in Figure Dl from LTR09. 



at the resonant radius tr. The radial width of the mode is 
Ar/ri? - 0.05. 

The linear growth of the Rossby wave is predicted to sat- 
urate when the azimuthal frequency of trapping of a fluid par- 
ticle in the trough of the wave ujt grows to a value equal to 
the growth rate Ui as discussed by LTR09. For the m — I mode 
the saturation level is estimated as \5p\/ p ^ lA{uji/ujr)'^{r/h) 
(LTR09). 
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